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FOREWORD 

The work described in this report is part of an alkali metal boiling and 
condensing heat transfer program conducted by the General Electric Company for 
the National Aeronautics and Space Administration under NASA Contract NAS 3- 
2528. The work was done under the technical management of Mr. F. E. Tippets, 
Missile and Space Division General Electric Company, and Miss Ruth N, Weltmann, 
Space Power Systems Division, NASA-Lewis Research Center. 
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ABSTRACT 

Forced convection heat transfer and pressure change data for condensation 

of potassium in vertical downflow inside single tubes are presented. These 

data, obtained at saturation temperatures from 1100°F to 1400°F, inlet vapor 

Mach numbers from about 0.1 to near 1.0 and local heat fluxes up to 
2 

300,000 Btu/hr-ft , provide a reasonable basis for thermal design of forced 
convection potassium condensers applicable to space power systems. 
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SUMMARY 


Tests to measure local condensing heat transfer coefficients and two- 

phase pressure change for potassium in forced convection were done in a 

stainless steel test rig at condensing temperatures from 1100°F to 140O°F, 

inlet vapor Mach numbers from about 0.1 to nearly 1.0, local vapor qualities 

of about 0.2 and 0.8, respectively, and local heat fluxes up to about 
2 

300,000 Btu/hr-ft . The tests were conducted using 36-inch long single 
tubes, cooled by sodium in counterflow, with the potassium in vertical 
downflow. Two different tube sizes were used, 3/8-inch and 5/8-inch ID. 

The tests were conducted both with and without inserts in the tubes. 

The local condensing potassium heat transfer coefficients were found to 
be relatively high, typically above 10,000 Btu/hr-ft 2 -°F. The data indicate 
that this value of the coefficient is a reasonably conservative design choice 
for forced convection condensing of potassium in tubes over the range of 
variables tested. 

The measured condensing heat transfer coefficients are generally lower 
than condensing potassium coefficients calculated by considering only the 
thermal resistance due to heat conduction through the liquid film, thus 
indicating the presence of an additional thermal resistance. It was possible 
to correlate the data by assuming the additional resistance was at the vapor- 
liquid interface and using an approach based on kinetic gas theory. This 
treatment resulted in an empirical value of 0.2 for the condensation coefficient 
associated with the resistance at the interface. 


Pressure change data was obtained by relating the measured potassium 
temperatures in the condenser to the corresponding saturation pressures. Although 
the pressure changes were quite small, in the order of less than 1-psi, reasonably 
good agreement was found between the experimental overall two-phase friction 
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pressure drop multipliers obtained from these measurements and those estimated 
using the conventional Martinelli flow model modified for potassium. 

The experimental results and associated analyses presented in this report 
provide a reasonable basis for thermal design of forced convection potassium 
condensers. 
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I INTRODUCTION 

Design of condensers for Rankine cycle space power systems using potassium 
as the working fluid requires knowledge of the heat transfer and fluid dynamic 
characteristics of condensing potassium in forced convection. At the time the 
work described in this report was initiated there were no data available on 
forced convection condensing of potassium in tubes. Therefore, the experi- 
mental program reported here was conducted to obtain basic data on potassium 
condensing in single tubes. 

The tests were done with the potassium in forced convection vertical 

downflow over the range of condensing temperatures from 1100°F to 1400°F, 

inlet vapor Mach numbers from about 0.1 to nearly 1.0 and local heat fluxes 

2 

up to about 300,000 Btu/hr-ft . The test results include pressure change 
data and measurements of local condensing heat transfer coefficients in 36-inch 
long tubes of two different sizes, 3/8-inch ID and 5/8-inch ID, both without 
inserts and with three different inserts in the 5/8-inch ID tube. 

Vertical downflow of the condensing potassium was chosen in order to 
eliminate the effect on the flow patterns in the condenser tube of having a 
component of the gravity force normal to the tube wall. The condenser tubes 
were cooled by sodium in counterflow outside the tube, in order to provide 
means for accurately measuring the total heat transfer rate and to have condensing 
heat flux distributions similar to those which would occur in a liquid metal 
cooled tube-in-shell condenser. Local condensing heat transfer coefficients were 
measured, rather than overall coefficients, in order to obtain better accuracy 
for the condensing coefficient determinations and to provide basic data which 
would be useful for development and design of space power condensers. 

The experimental results and associated analyses are presented in detail 
in the following sections of this report. 
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II SUMMARY OF LITERATURE SURVEY 

In general, condensation heat transfer can be considered to occur in one 
or a combination of three different modes, as follows: 

(1) Droplet formation in the vapor stream 

(2) Dropwise condensation on the cooling surface 

(3) Film condensation on the cooling surface 

Due to its wetting characteristics, potassium probably condenses as a 
liquid on the cooling wall in forced convection flow. Therefore, the following 
discussion of the literature will treat only this mode of condensation. 


Until recent times, data on condensation heat transfer has been available 
only for fluids with high Prandtl numbers. For these fluids the analytical 
prediction of Nusselt for a falling film in a stagnant vapor has been found 
to generally apply (Reference 12). However, more recently available data on 
condensing heat transfer for liquid metals of low Prandtl number indicate that 
the simple liquid film analysis proposed by Nusselt is not adequate for these 
fluids (References 13-20). 

A number of analytical investigations (References 21-29) have been undertaken 
to modify Nusselt’ s model to include such effects as vapor shear, film turbulence 
and fluid property variation across the liquid film. However, none of these 
treatments, for which it is assumed that the liquid film is the only thermal 
resistance, succeed in bringing the analytical predictions into agreement with 
measured condensing heat transfer coefficients for the low Prandtl number fluids. 

Some evidence of the existence of a substantial additional thermal resistance 
besides that associated with the liquid film has been presented by Rohsenow and 
Sukhatme for mercury condensing on the outside of a 3/4-inch OD x 6-inch long tube 
in a stagnant vapor chamber (Reference 16). Their results for mercury condensing 
on the outside of this vertical tube were: 
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(1) The measured liquid film thicknesses are in reasonable agreement 
with those predicted by the Nusselt model . 

(2) The measured condensing heat transfer coefficients are as much 
as two orders of magnitude lower than those predicted by the 
Nusselt theory. 

On the basis of these observations, Sukhatme and Rohsenow concluded that 
for mercury and other low Prandtl number fluids there is a significant vapor 
phase thermal resistance, which is additive to the liquid film resistance. 

This vapor phase thermal resistance also may be present with high Prandtl 
number fluids, but for these fluids it would be negligible in comparison with 
the thermal resistance of the liquid film. 

As pointed out by Wilhelm (Reference 30) , the presence of a non- 
condensable gas can lower the condensing heat transfer coefficient. Thus, this 
is another possible cause of the large difference between measured condensing 
heat transfer coefficients for low Prandtl number fluids and those calculated 
using the liquid film analyses. 
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Ill TEST APPARATUS 


Facility 

An isometric drawing of the test facility is presented in Figure 1. 

Figure 2 shows the facility during initial stages of construction and 
Figure 3 shows the completed facility prior to initial operation. The 
entire facility, except the test section, is constructed of Type 316 
stainless steel, which was selected on the basis of its strength and 
compatibility with alkali metals at temperatures to 1600°F. 

Potassium vapor was generated in a pot boiler heated by electrical 
immersion heaters. The pot boiler, shown in Figure 4, was constructed of 
a 12-inch Schedule 40 pipe tee and caps, with a flat plate serving as the 
flange for the heaters. Twenty immersion heaters, each rated at 4.7 KW 
were installed in the boiler, as shown in Figure 5. These heaters were 
electrically connected so that no single heater could exceed a 3.6 KW 
electrical input. Two types of artificial nucleating sites were employed 
in the boiler. The first type, referred to as a "hot finger", consisted 
of a 1-inch OD Type 316 stainless steel bar containing a 0.040-inch OD hole. 

Three of these "hot fingers" were located in the bottom of the boiler. 

Independent electrical heaters maintained the "hot fingers" above the 
potassium saturation temperature. A second type of artificial nucleator, which 
was made of one-inch long sleeves shrunk onto the bottom row of immersion heaters, 
was also provided. No measurable effect of either of these types of nucleators 
on loop operation was noted. 

As indicated in Figure 1, the vapor generated in the boiler flowed through 
an 8-foot vertical length of two-inch Schedule 40 pipe, a 1^-inch throttle valve, 
a 10-foot horizontal run of 1-inch Schedule 40 pipe, and then passed down through 
the sodium-cooled test section where it was condensed. During test operation, the 
liquid-vapor interface was maintained below the test section in the head tank. 
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After exiting from the test section and head tank, the liquid condensate 
passed through a thermal calorimeter, an induction electromagnetic pump 
and an electromagnetic flowmeter into the boiler. All pipe lines were 
wrapped with Inconel sheathed Chromel "a” heating wire for preheating. 

The sodium loop, an all-liquid loop, was used to cool the test condenser. 
The sodium flowed vertically upward through the test section annulus where it 
was heated by the condensing potassium, flowed through 8-foot lengths of 1-inch 
and 2-inch Schedule 40 pipe, respectively, and then passed through a conduction 
electromagnetic pump and an electromagnetic flowmeter. The sodium then passed 
through a tube-in-shell air cooler where heat was rejected, after which the 
sodium returned to the test section. A dump tank was provided in the loop for 
storing the sodium during shutdown. 

For purification of the potassium, a bypass hot trap, shown in Figure 6, 
was installed in parallel with the test condenser. The hot trap, which 
contained approximately 4 lbs of zirconium gettering material, provided hot 
trapping to remove oxide from the potassium during flushing of the loop prior 
to the beginning of operation. In the sodium loop, a diffusion cold-trap was 
provided. This cold-trap was in the form of a cold leg between the sodium pump 
and dump tank to allow oxide precipitation during operation. 

Test Section 

A schematic drawing of the condensing test section used in this series of 
experiments is shown in Figure 7. Potassium condensation occurred inside a 
thick walled nickel tube which was cooled by a countercurrent sodium flow in 
the surrounding annulus. Operation was with the condenser tube axis vertical 
and the potassium in vertical downflow. 

During this investigation, experimental data were obtained for five 
condensing geometries, which are listed below in the chronological order of 
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testing: 


Geomet ry 

5/8-inch ID tube with no insert 

5/8-inch ID tube with an uni st rument ed tapered pin insert 
3/8-inch ID tube with no insert 

5/8-inch ID tube with an instrumented helical insert 
5/8-inch ID tube with an instrumented 1/4-inch OD tubular insert 

The three inserts which were tested are shown in Figure 8. 

The condenser tubes were thick-walled and were made of high-purity nickel 
(99.95%+) which was selected for its high thermal conductivity in combination 
with its compat ability with the alkali metal test environment. The tubes 
had an active heat transfer length of 36-inches (Figure 7). The potassium flow 

passage was concentric with the outside surface diameter of the tube to within 
0. 003-inch. 

To calculate the condensing heat transfer coef f icient s, measurements of 
the local heat flux, inner wall temperature and vapor temperature are required. 
These measurements were made at two stations in the condenser, one at about 20% 
and the other at about 80% of the condensing length measured from the potassium 
vapor inlet (Figure 7). To determine the inner wall temperature and heat flux 
at these two measuring stations, thermocouples were located in holes drilled 
by an electrical discharge process from each end of the thick-walled nickel 
tube parallel to the axis of the tube. There were five of these holes at 
each end of the nickel tube, located as shown in Figure 9. 

The drilled thermocouple holes were 0.050-inch in diameter and 9-inches 
deep, and their radial location was held within a tolerance of 0.004-inch over 
their length, as determined by radiographic examination. Two of the five holes 


Number 

1 

2 

3 

4 

5 
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at each measuring station were located close to the tube inside surface, one 
hole was located at an intermediate radial position, and two holes were located 
close to the tube outside surface (Figure 9). One thermocouple was placed in 
each of the thermocouple holes at the depths for each measuring station indicated 
in Figure 7. The axial positions indicated in Figure 7 are estimated to be 
correct within less than + 1/16-inch error. 

Potassium vapor temperatures were measured only at the test section inlet 
and outlet with the first three test section geometries. However, with test 
section geometries No. 4 and No. 5, thermocouples were located in the test 
section inserts along the condensing length to measure directly the potassium 
axial temperature distribution. 

In addition to the thermocouples located within the nickel tube wall and 
in the potassium stream, thermocouple wells were provided to measure the sodium 
inlet and outlet bulk temperature. Figure 10 shows the sodium inlet end of the 
test section with three thermocouple wells spaced 120° apart for measuring the 
inlet bulk sodium temperature. A similar arrangement was provided at the sodium 
exit . 


Figure 11 shows one end of the 3/8-inch inside diameter nickel condenser 
tube with a stainless steel adaptor pipe brazed into place for connection of 
the nickel condenser tube to the potassium loop piping. Arrangements identical 
to this were used at both ends of each test section. In addition to the adaptor 
pipe arrangement shown in Figure 11, each end of the nickel condenser tube was 
fitted with a 2-inch OD stainless steel tube which was welded to the outside 
surface of the nickel tube at each end and was used to join the nickel tube 
to the surrounding sodium-containing shell. The space between the adaptor pipe 
and the 2-inch OD tube connecting the condenser tube to the shell formed an 
annulus which was open to the atmosphere, through which the tube wall thermo 
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couples were Installed. This arrangement; which can be seen in Figure 7; 
avoided crossing either the sodium or potassium streams with the wall thermo- 
couple leads; thereby avoiding mechanical complexity and undesirable flow 
disturbances. 

The test section outer shell was constructed of 2^-inch Schedule 80 pipe 
with two enlarged end fittings made from standard pipe reducers and flat plates. 
Each end of the test section was internally baffled on the sodium side with a 
perforated disc. Concentricity between the test section shell and the thick 
walled nickel condenser tube was maintained by small spacers which were 
attached to each end of the nickel tube. Flexibility to accommodate lengthwise 
differential thermal expansion between the stainless steel shell and the nickel 
tube was provided by a guided 3^-inch OD bellows located at the bottom end of the 
outer shell. 

Figure 12 shows a condenser tube with the adaptor pipes attached and shows 
an assembled condensing test section. The condenser tube is shown adjacent to 
the assembled test section to indicate its relative length and location within 
the test section shell. 

Figure 13 shows the test section installed in the facility before covering 
it with thermal insulation. As shown; in the Figure; the test section was 
surrounded by a structural framework to minimize bending and to maintain a 
vertical orientation. The supporting structure allowed the test section to 
expand vertically but prevented deviation from vertical orientation of more than 
+ 1/16-inch in the 28-inch length between the supports. The shell thermocouples 
shown in Figure 13 were used during initial operation with Test Section No. 1 to 
check on the axial symmetry of the sodium flow. For subsequent operation; two 
thermocouples were attached to the shell at the sodium inlet; at the sodium outlet 
and in the middle of the test section. 
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The first of the three inserts tested consisted of a non- instrumented 
tapered pin which had diameters of 1/32-inch and 3/8-inch at the test section 
inlet and outlet, respectively (Figure 8-a) . This insert was used to simulate 
a tapered tube flow geometry by providing a gradually reducing flow cross- 
section to the condensing potassium stream. The insert was held concentric with 
the inside diameter of the condensing tube by wire prongs welded at two positions, 
approximately 11-inches and 30-inches, respectively, from the test section inlet. 

The tapered pin was also secured in the axial direction by a hairpin shaped yoke 
attached to the upstream end of the insert and tack-welded to the pipe inside 
surface above the test section inlet. 

The second insert used was an instrumented helical insert of pitch-to- 
diameter ratio equal to 6 (Figure 8-b) . This insert consisted of a single 
helical ribbon wound around a 1/4-inch outside diameter centerbody tube 
(Figure 14). The centerbody tube penetrated the potassium loop piping through 
a fitting located above the test section inlet, as shown in Figure 14. Twelve 
thermocouples were located in the insert centerbody tube. Three of the thermo- 
couples were located upstream of the active condensing heat transfer length, 
and the other nine thermocouples were placed at positions about 5-inches apart 
along the active condensing length (two thermocouples were located together at 
the potassium outlet). The loop piping penetration for the centerbody tube was 
made at an enlargement in the potassium loop piping, which consisted of two 
eccentric reducers welded together back to back. This arrangement was used for 
two purposes; first, to permit a gradual long radius bend in the insert tube 
which made insertion and removal of the thermocouples possible without requiring 
bending of the thermocouples, and secondly, to provide a penetration location 
which had the proper angle and accessibility for satisfactory welding. 

The third insert was similar in arrangement to the one shown in Figure 14, 
except that instead of being a helix it was a straight 1/4-inch OD tube 
(Figure 8-c) . This straight tubular insert was held concentric with the test 
section by means of wire prongs and guide vanes. The tubular insert also contained 
twelve internal thermocouples, which were distributed in a manner similar to that 
used for the helical insert, to measure the axial temperature distribution of the 
condensing potassium. 
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Figure 1, Flow Diagram of Condensing Test Facility 
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Figure 3. Condensing Test Facility Prior to Initial Startup 
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Figure 4. Photograph of Pot Boiler Showing Arrangement of 
Immersion Heaters and Artificial Nucleators 
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Figure 9* End View of 3/8-inch ID Condenser Tube Showing 
Drilled Thermocouple Holes 
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Figure 10. End View of Test Condenser Shell Showing Sodium Inlet 
Thermocouple Wells 







Figure 11. End View of 3/8-inch ID Condenser Tube Showing Potassium 
Loop Adaptor Pipe and Wall Thermocouple Holes 
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Sodium Outlet 


Test Section 


Figure 13 


Test Condenser Installed In Condensing Test Facility 
(Shell thermocouples shown were used during initial 
operation with Test Section No. 1 to check symmetry 
of sodium flow) . 




- 26 - 


IV EXPERIMENTAL TECHNIQUE 


Test Procedure 

During two-phase operation, the independent parameters in the 
Condensing Test Facility were: 


1. Boiler power, q fe 

2. Air flow rate to sodium air cooler, 

3. Sodium flow rate, W^ a 

4. Auxiliary heat input 

The dependent parameters are: 

1. Potassium flow rate, W R 

2. Potassium vapor temperature, T R 

3. Sodium temperature, T^ a 

4. Test section inlet quality, 

Auxiliary heat input and heat losses are assumed to be negligible, and 
the potassium is assumed to be 100% saturated vapor at the test section inlet. 
With these assumptions the potassium vapor flow rate is related to the boiler 
power, as expressed by Equation (1) and Sketch a. 
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Neglecting heat losses, the boiler heat input is equal to the heat 
rejected to the sodium coolant in the test condenser, as expressed by 
Equation (2). 


q C W Na Cp Na (T NaO T NaI 5 


( 2 ) 


The coupling between the sodium and potassium loops occurred in the 

test condenser and can be expressed in terms of U the overall heat transfer 

coefficient, A the test condenser heat transfer area, (T. - T ) the log 
c K Na LM 

mean temperature difference, and q the test section heat input, as given 
by Equation (3). 


= U A r 


(T - T ) 

K Na'LM 


(3) 


The heat transfer resistance of the thick-walled nickel tube was large 
compared to the sodium and potassium thermal resistances. Thus, the overall 
heat transfer coefficient was approximately constant during testing. The 
liquid level was maintained below the test section exit. Thus, the condensing 
heat transfer area was constant and was equal to the total area of the test 
condenser. Therefore, the coupling between the sodium and potassium loops 
expressed by Equation (3), was as illustrated by Sketch-b and Sketch-c. 



Sketch-b 
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The heat gained by the sodium in the test section was rejected to air 
in the sodium loop air cooler (Figure 1). The heat rejection rate is given 
by Equation (4), for which it is assumed the air-side heat transfer coefficient 
h A 1S controlling. Sketch-d illustrates the functional relationship between the 
heat rejection rate and the sodium temperature. 

q C = ^ = h A A AC (T Na ' T A>LM (4) 



Sketch-d 
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The sodium flow rate was maintained constant at its maximum value 
during testing. Thus, the sodium flow rate was not employed as a test 
parameter. 

The system response to an increase of boiler power will now be 
discussed. An increase of boiler power results in a corresponding increase 
in potassium vapor temperature, in order to enable rejection of the additional 
heat input in the test section. With a constant air flow rate to the sodium 
cooler, the sodium temperature also increases to allow rejection of the additional 
heat rate to the air, as indicated in Sketch-d. When steady-state is reached, 
the potassium temperature will have increased more than the sodium temperature 
increase, since, as indicated in Sketch-b the test section log mean temperature 
difference between the potassium and the sodium must also increase in order to 
allow an increase in the condenser heat transfer rate. The potassium flow rate 
will have increased also, due to the increased boiler power, as shown in 
Sketch-a. If the air flow rate is increased at fixed boiler power, the sodium 
temperature will decrease, thereby lowering the potassium condensing temperature, 
as indicated in Sketch-c. Thus, at fixed boiler power the potassium condensing 
temperature can be controlled by adjusting the air flow rate to the sodium 
cooler. These characteristics of the system formed the basis of the condensing 
test procedure. The sequence of steps used in the operating procedure was as 
follows: 


1. The potassium loop was first evacuated at 600°F to a pressure of 
approximately 25 microns. The difference between liquid and vapor temperatures 
in the boiler was used to check for the partial pressure of inert gases in the 
system. If a substantial difference (greater than 10°F) between those 
temperatures was obtained, the loop was re-evacuated and the process was repeated 
to the extent required. 
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2. Next, sodium flow was established, a boiler power was set, and 
the temperature of the sodium was adjusted using the air cooler and line 
heat to bring the potassium vapor at the test section inlet to the desired 
temperature. 

3. Tests were conducted at constant boiler power by systematically 
varying the sodium temperature to obtain data for potassium vapor inlet 
temperatures over the range from 1100 to 1400°F. 

4. When a group of tests was completed at a given boiler power, the 
boiler power was changed, and then step-3 was repeated. 

This procedure was repeated to provide condensing heat transfer data, 
with condenser heat transfer rate and condensing temperature as parameters, 
for each of the five test section geometries. Figure 15 shows the range of 
operating capability of the Condensing Test Facility, within which the 
condensing tests were conducted. Line CD is the 1600°F temperature limit of 
the facility and line AE sets the lowest temperature of interest at 1100°F. 

Line AB sets the boiler power limit for choked flow in the test condenser. 

Line BC sets the maximum boiler power at 50 KW. Line DE sets the minimum 

net boiler power at 4 KW for a minimum condensing heat flux of 30,000 Btu/hr-ft 2 

for a 36-inch long 5/8-inch ID condenser. 

Prior to the beginning of condensing test operation for Test Set No. 1, 
liquid tests were made to obtain operating experience with the test facility 
and to gain familiarity with the test section instrumentation. Potassium liquid 
heat transfer coefficients were calculated from these data. These data are 
presented in Appendix A. 

Flow Rate Measurement 

The potassium and sodium flow rates were measured using electromagnetic 
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flowmeters. A calorimetric method was used to calibrate the flowmeters, as 
described in Appendix B. 

Temperature Measurement 

Sheathed, capped, 0.040-inch diameter, chrome 1-alumel thermocouples were 
selected for the temperature measurements, on the basis of their stability and 
relatively high emf output in the 1100-1400°F temperature range. The high 
output of the chromel-alumel thermocouples, approximately 23 micro-volt s/°F, 
minimized temperature errors due to instrumentation noise. The instrumentation 
noise was generally less than + 3 micro-volts, which correspond to + 0. 13°F. 

During the early stages of the program it became evident that the 

condensing heat transfer coefficient for potassium was relatively high, in the 

2 

order of 10,000 Btu/hr-ft -°F. Thus, accurate calibration of the test section 
thermocouples became necessary. For the calculation of the condensing heat 
transfer coefficient a difference between the vapor and wall temperatures is 
needed. Therefore, it was not necessary to calibrate each thermocouple on 
an absolute basis, but only to obtain accurate in-place calibration of the test 
section thermocouples relative to each other. 

Both before and after each test set, the thermocouples in the potassium 
fluid and in the nickel tube wall were calibrated in-place relative to one of 
the potassium thermocouples, at temperatures of 1200 and 1300°F, as discussed 
in detail in Appendix C. On the basis of these calibrations each thermocouple 
was assigned a temperature correction as a function of temperature. The sodium 
well thermocouples were also calibrated in-place relative to one of the inlet 
sodium well thermocouples using a technique described in Appendix C. 

Early calibration attempts were hampered by inhomogeneities in the 
chromel-alumel sheathed thermocouples, which were caused by small bends made 
during installation. These inhomogeneities were detected by moving a soldering 
iron at about 700°F along the exposed thermocouple length, after installation, 
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and observing the thermocouple reading. In the early stages of the program 
it was not uncommon to obtain outputs due to inhomogeneities as great as 
500 micro-volts (20°F) . Thermocouples for which the homogeneity check 
indicated an output greater than 50 micro-volts were replaced. In addition, 
individual ice junctions were used for each thermocouple to reduce errors 
caused by junction emf's. 

As a further effort to assess temperature measurement accuracy, the 
magnitude of distortion of the temperature field in the thick-walled nickel 
tube caused by the axial thermocouple holes was estimated analytically. The 
results of this study are discussed in Appendix D. 

To calculate the condensing heat transfer coefficient, the potassium 
saturation temperature, heat flux, and inner wall temperature had to be 
determined. During reduction of the condensing data, the inner wall temperature 
and heat flux at each measuring station were obtained from the nickel tube 
wall radial temperature profile. The method used is discussed in Appendix E, 
which describes the condensing data reduction in detail. Figure 16 shows two 
typical radial temperature profiles (obtained during Test Set No. 5). As 
discussed in Appendix E, a least squares procedure was employed to fit the 
measured temperature within the nickel tube wall using the integrated Fourier 
conduction equation. The local inner wall heat flux, q/A ± , was obtained from 
the slope of the radial wall temperature profile. 

For the five test sets, two different methods were used to determine 
the local potassium temperature at the measuring stations. For the first 
three test sets, which included the 5/8-inch and 3/8-inch ID plain tubes 
without inserts and the 5/8-inch ID tube with tapered pin insert, only the 
potassium inlet and outlet temperatures were measured. For these data the 
local potassium temperature at the measuring stations was estimated by linear 
interpolation between the measured inlet and outlet temperatures, with correction 
for friction pressure drop in the inlet pipe from the temperature measuring 
point to the test section inlet. For the last two test sets, which consisted of the 
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5/8-inch ID tube with the instrumented helical insert (Test Section No. 4) 
and with the instrumented tubular insert (Test Section No. 5), the local 
potassium temperature was measured directly with the insert thermocouples. 

In Figure 17 are shown two potassium axial temperature distributions 
and corresponding saturation pressures at an inlet temperature of 1150°F 
for two mass velocities which were measured in the 5/8-inch ID tube with 
tubular insert. The pressure profiles indicate that in both cases the 
frictional pressure drop near the test section inlet exceeds the pressure 
rise due to momentum reduction, but is less than the momentum pressure rise 
near the test section exit. A net pressure increase is indicated over the 
condensing length for both cases.. A larger pressure rise across the test 
section was obtained at the higher mass velocity. 

In Figure 18 are shown two potassium axial temperature distributions 
and corresponding saturation pressures for inlet temperatures of 1300°F 

2 

and 1400° F at a superficial mass velocity of approximately 13 lb/sec-ft , 
which were also measured in the 5/8-inch ID tube with tubular insert. For 
the lower inlet temperature, a larger momentum pressure increase was obtained. 
This is caused by the larger inlet vapor velocity at the lower temperature 
for the same flow rate. The precision of the insert thermocouples is indicated 
by examination of Figure 18b. As can be seen, the total temperature variation 
between the test section inlet and outlet was less than 2°F, and the maximum 
scatter of the measured temperatures is less than 0.5°F at 1400°F. Included 
in Figure 18 are the measurements of the last two of the three thermocouples 
which were located in the vapor inlet pipe upstream of the condensing heat 
transfer section. These thermocouples show that the pressure in the inlet 
vapor pipe decreases as a function of length, due to friction pressure drop. 

Figure 19 shows two typical potassium temperature distributions measured 
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using the instrumented helical insert (Test Set No. 4), one for a mass 
velocity of 7.4 lbs/sec-ft 2 and the other for a mass velocity of 14 lbs/sec-ft 2 
at an inlet vapor temperature of 1200°F. As can be seen from Figure 19, there 
was a net pressure drop across the condenser with the helical insert, due to 
the friction pressure drop being larger than the momentum pressure rise. 
Increasing the mass velocity tended to increase the net pressure drop 
(compare Figure 19b with Figure 19a). 
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Figure 17. Potassium Axial Temperature Distributions and 
Corresponding Saturation Pressures For £wo 
Different Mass Velocities At About 1150 F Inlet 
Temperature (5/8-inch ID tube with l/4-inch 0D 
Instrumented Tubular Insert) 
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Figure 18. Potassium Axial Temperature Distributions and 

Corresponding Saturation Pressures For Two P 

Different Inlet Temperatures at G = 15 lb/sec-fr 
(5/8-inch ID Tube with l/4-inch OD Instrumented 
Tubular Insert) 
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Figure 19- Potassium Axial Temperature Distributions and 
Corresponding Saturation Pressures for Jwo 
Different Mass Velocities at About 1200 F Inlet 
Temperature (9/8-inch ID Tube with Instrumented 
Helical Insert, p/D. = 6) 





V CONDENSING POTASSIUM HEAT TRANSFER RESULTS 


A summary of the test conditions used with each of the five condenser 
test section geometries is given in Table 1. The experimental condensing 
heat transfer data obtained are given in Tables 2 to 6. The data reduction 
procedures employed are described in Appendix E. A discussion of estimated 
errors in determination of the condensing heat transfer coefficients is given 
in Appendix F. 

For purposes of discussion, the results from Test Set No. 4 (5/8-inch ID 
tube with instrumented helical insert, P/D. = 6) and Test Set No. 5 (5/8-inch ID 
tube with instrumented 1/4-inch OD tubular insert) are presented together. These 
data, given in Tables 2 and 3, are the most reliable of the data obtained for 
the five test sections because the local potassium temperature was measured 
directly using thermocouples in the inserts. The results from Test Sets No. 1, 

2 and 3, for which the local potassium temperatures were estimated by linear 
interpolation from measurements at the test section inlet and outlet, are 
presented as a second group. 

Results From Test Sets No. 4 and 5 

The heat transfer results obtained using an instrumented helical insert in 
a 5/8-inch ID tube (Test Set No. 4) and an instrumented 1/4-inch OD tubular 
insert in the same 5/8-inch ID tube (Test Set No. 5) are presented in Figures 20 
and 21 and in Tables 2 and 3. The two sets of data are presented together for 
comparison in Figure 22. Calculated values of the condensing heat transfer 
coefficients at 1300°F, calculated using Nusselt's laminar film model (Reference 12) 
and Seban's turbulant film model (Reference 21), are also shown in the Figures 
for comparison with the data. 

The tubular insert data from Test Set No. 5 are based directly on the local 
measurements of the potassium temperature obtained from the insert thermocouples. 
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The helical insert data from Test Set No. 4 are based on the local measurements 
of the potassium temperature obtained from the insert thermocouples, corrected 
to account for the increased pressure and saturation temperature at the tube 
wall due to the radial acceleration caused by the helical insert. The analysis 
on which this correction is based is given in Appendix E. 

The local vapor qualities for each of the data points taken with an 
instrumented tubular insert in Test Set No. 5, listed in Table 3, were calculated 
by energy balance using the measured fluid temperatures and flow rates, as 
discussed in Appendix E. The approximate local vapor qualities at the top 
and bottom measuring stations, respectively, indicated in Table 2 for the helical 
insert data of Test Set No. 4 were estimated by linear interpolation assuming 
100% vapor at the test section inlet and 0% vapor at the test section outlet. As 
discussed in Appendix E, this procedure results in an estimated error of about 
5% in local quality at the measuring stations, compared to vapor qualities 
calculated by energy balance. This amount of error is not significant in 
determination of the condensing coefficients. A more detailed calculation of 
the local vapor quality is not justified for the helical insert data because of 
the uncertainties inherent in correcting for radial acceleration effects due to 
the helical flow. 

Examination of Figures 20-22 shows that: 

(1) The condensing heat transfer coefficients are in the order of magnitude 

o 

of 10,000 Btu/hr-f t -°F or more (Figures 20 and 21). 

(2) The data from the bottom measuring stations {t / = 48) where the 
vapor quality is least (X ~ 0.20) have slightly lower heat transfer coefficients 
than those from the top measuring station (/^/D^ = 10) where the vapor quality is 
higher (X ,=y 0.88). As can be seen from the Figures, this experimental trend is in 
general agreement with the trend of reduced heat transfer coefficient with increased 
liquid film Reynolds Number predicted by theory (References 12 and 21). 

(3) The measured condensing heat transfer coefficients are less than those 

calculated by the Nusselt and Seban models. 
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(4) Agreement between the helical insert data from Test Set No. 4 and 
the tubular insert data from Test Set No. 5 is within the scatter of the dat 
(Figure 22). 


The local condensing heat transfer coefficient data taken with an in- 
strumented 1/4-inch OD tubular insert in a 5/8-inch ID tube (Test Set No. 5) 
are plotted in Figure 23 versus local potassium saturation temperature, mass 
velocity and local heat flux. As can be seen in Figure 23, the condensing 
heat transfer coefficient is approximately independent of heat flux, mass 
velocity and saturation temperature over the ranges of conditions tested. 


Three of the data points shown in Figure 23 were taken at condensing 
temperatures in the range from about 880°F to 1060°F. As shown in Figure 23, 
the heat transfer coefficients for these three data points are significantly 
less than for the other data, which were taken at saturation temperatures in 
the range from 1100°F to 1400°F. It is possible that this is due to a trend 
of decreased heat transfer coefficient with decreasing temperature which might 
exist at the lower temperatures but is not evident in the higher range of 
temperatures in which the rest of the data were taken. It is also possible 
that these three points have lower coefficients because of some undetermined 
experimental error which occurred when they were being taken. Excepting these 
three lower temperature data points, nearly all the other Test Set No. 5 data 
have condensing coefficients that are larger than 10,000 Btu/hr-ft 2 -°F, as 
shown in Figure 23. 


The characteristic of the measured condensing heat transfer coefficients 
being generally less than those calculated using the liquid film models of 
Nusselt (Reference 12) and Seban (Reference 21), as shown in Figures 20-22, 
suggests the possibility of an additional thermal resistance besides that 
associated with conduction through the liquid film. An investigation of this 
was done by treating the Test Set No. 5 data (5/8-inch ID tube with instrumented 
1/4-inch OD tubular insert) using an analytical model based on the concept of 
there being a vapor phase thermal resistance that is additive to the resistance 
due to conduction through the liquid film, as suggested by Rohsenow and Sukhatme 
(Reference 16). The analysis is presented in detail in Appendix G. Major steps 
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in the treatment of the Test Set No. 5 data to determine the vapor phase 
thermal resistance and corresponding vapor phase coefficients are as follows. 

First, the liquid film thickness 6 was calculated using Dukler' s model 
(Reference 28) for each of the Test Set No. 5 data points. Using these 
values of 6 the heat transfer coefficient corresponding to conduction through 
the liquid film was calculated employing Equation (5). 


h f = k f //6 

The vapor-phase thermal resistances (1/10 were then determined by 
subtracting the calculated liquid film thermal resistance (1/h^) from the 
experimental values of the total condensing thermal resistance (l/h c ) measured 
for each of the Test Set No. 5 data points, using Equation (6). 



In order to compare with theory the values of h^ calculated from the 
experimental heat transfer data using Equation (6), an analysis was done 
based on kinetic gas theory following a procedure suggested by the work of 
Schrage (Reference 32). The analysis, given in Appendix G, results in Equation (7) 
for the vapor phase coefficient h v as a function of vapor saturation properties, 
absolute temperature and condensation coefficient O . 


h = 
v 


a X 

c 


M gp 
2 n R T 


- £-) 
aT 2T 


(7) 


The vapor-phase heat transfer coefficients obtained from the Test Set No. 5 
data using Equation (6) are presented in Figure 24, together with calculated 
values from Equation (7) for 0^ = 0.2. The selected value of ° c = 0.2 gives an 
approximate best fit to the data (average deviation less than + 25%). The 
vapor-phase heat transfer coefficients obtained for each of the Test Set No. 5 
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each data point 


data points and the corresponding values of or calculated for 

c 

using Equation (7) are listed in Table 4. 

As shown in Figure 24, the selected value of ^ = 0 .2 results in fair 
agreement between the vapor-phase heat transfer coefficients calculated using 
Equation (7) and those obtained from the Test Set No. 5 data using Equation (6). 
Vapor-phase coefficients obtained from the data of Engelbrecht (Reference 19) 
using Equation (6), shown in Figure 24, are in order-of-magnitude agreement with 
values calculated using Equation (7) with ^ = 0. 2 for saturation temperatures 
down to 700°F. A trend of decreasing vapor-phase heat transfer coefficient with 
decreased saturation temperature is calculated using Equation (7) with <7 =0.2, 

as shown in Figure 24. 

Results From Test Sets No. 1, 2 and 3 

The heat transfer results obtained using a 5/8-inch ID tube without insert 
(Test Set No. 1), a 5/8-inch ID tube with non-inst rument ed tapered pin insert 
(Test Set No. 2) and a 3/8-inch ID tube without insert (Test Set No. 3) are 
presented in Figures 25-27 and Tables 5-7, respectively. Calculated values of 
the condensing heat transfer coefficients at 1300°F, calculated using Nusselt's 
laminar film model (Reference 12) and Seban 1 s turbulent film model (Reference 21), 
are also shown in the Figures for comparison with the data. 

For Test Sets No. 1, 2 and 3 the local potassium temperatures at the 
measuring stations, on which the local condensing heat transfer coefficients are 
based, were estimated by linear interpolation between the measured potassium inlet 
and outlet temperatures, with correction for pressure drop in the vapor inlet pipe 
between the upstream bulk thermocouple station and the condensing section inlet. 

As discussed further in Appendix H, this procedure for estimating the local 
potassium temperature at the measuring stations probably contributed to the 
increased scatter evident in the data obtained from Test Sets No. 1, 2 and 3. 

The approximate local vapor qualities at the measuring stations indicated 
in Tables 5-7 for Test Sets No. 1, 2 and 3, respectively, were estimated by 
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linear interpolation assuming 100% vapor at the test section inlet and 0% 
vapor at the test section outlet. Based on the calculated quality distribution 
comparisons given in Appendix E this procedure results in an estimated error 
of about 5% in local quality at the measuring stations. This amount of error 
is not significant in determination of the condensing coefficients. 

Examination of Figure 25 shows that the heat transfer coefficient data 

from Test Set No. 1 (5/8-inch ID tube without insert) have more scatter but are 

in order-of-magnitude agreement with the Test Set No. 5 data (5/8-inch ID tube 

with instrumented 1/4-inch OD tubular insert) shown in Figure 21. Most of the 

2 

Test Set No. 1 data have coefficients larger than 10,000 Btu/hr-ft -°F, but all 
of the coefficients are lower than values calculated using the Nusselt laminar 
film model (Reference 12) and the Seban turbulent film model (Reference 21), as 
shown in Figure 25. 

The Group-D data listed under Test Set No. 1 in Table 1, comprised of 
forty-one data points taken in a 5/8-inch ID tube without insert, are not 
included in Figure 25 and Table 5. This group of data was taken at approximately 
the same conditions as the other data from Test Set No. 1 given in Figure 25 and 
Table 5 and is reported in Reference 8. The majority of the Group-D data had 
indicated coefficients that were negative due to excessive thermocouple errors 
and possible mal-distribution of the sodium flow in the condenser shell, as 
discussed in Reference 8. Thus, the Group-D data are not useable. These sources 
of experimental error were corrected for the other tests. 


The data from Test Set No. 2 (5/8-inch ID tube with non-inst rumented tapered 
pin insert), are shown in Figure 26. The coefficients have more scatter and are 
higher than the data from Test Sets No. 1 and 5, and they are higher than values 
calculated using the Nusselt model and the Seban model. Examination of Table 6 
shows that about half of the data points taken during Test Set No. 2 have indicated 
coefficients that are negative, primarily at the top measuring station (not plotted 
in Figure 26) , indicating excessive experimental errors in the data. Thus, the 
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Test Set No. 2 data are not reliable and in consequence no definite conclusions 
can be drawn regarding the effect of the tapered pin insert on the condensing 
heat transfer coefficient. 

Figure 27 shows the data from Test Set No. 3 (3/8-inch ID tube without 
insert). The coefficients have large scatter, both above and below values 
calculated using the Nusselt and Seban models, especially at the top measuring 
station ( ,2/Ih = 21). Examination of Table 7 shows that some of the data from 
the top measuring station have indicated coefficients that are negative (not 
included in Figure 27). Due to the large scatter in the data shown in Figure 27, 
no definite conclusions can be drawn regarding the effect on the condensing 
coefficient of a smaller tube size (3/8-inch ID) in comparison with coefficients 
measured in the 5/8-inch ID tubes (Test Sets No. 1 and 5). 
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Condensing Heat Transfer Coefficient, h , Btu/hr-ft 



Figure 20. Local Condensing Heat Transfer Coefficients From Test 
Set No. 4, 5 / 8 -inch ID Tube with Instrumented Helical 
Insert (p/D. = 6 ) 
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Figure 21. Local Condensing Heat Transfer Coefficients From 

Test Set No. 5, ^/Q-lnch ID Tube with Instrumented 
l/4-inch 0D Tubular Insert. 
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Figure 22. Local Condensing Heat Transfer Results from 
Test Sets No. 4 and 5/8-inch ID Tube with 
Instrumented Inserts 
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Figure 23. Variation of Local Potassium Condensing Heat Transfer 
Coefficient with Saturation Temperature, Mass Velocity 
and Heat Flux (5/8-inch ID Tube with l/ 4 -inch OD 
Instrumented Tubular Insert) . 
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Figure 25. Local Condensing Heat Transfer Coefficients from 
Test Set No. 1 , 5/8-inch ID Tube Without Insert 
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Figure 26. 


Local Condensing Heat Transfer Coefficients From 
Test Set No. 2, 5 / 8 -inch ID Tube With Non-Instrumented 

Tapered Pin Insert. 



Condensing Heat Transfer Coefficient, h , Btu/hr-ft 



Figure 27. Local Condensing Heat Transfer Coefficients From 
Test Set No. 5, 5 / 8 -inch ID Tube Without Insert. 







VI CONDENSING POTASSIUM PRESSURE CHANGE RESULTS 


In this section the experimental two-phase pressure change data 
obtained during the condensing tests are presented and briefly compared 
with predictions using the homogeneous flow model (K = 1) and the 
Martinelli model (Reference 31) modified for application to potassium 
(Reference 6). The experimental data presented include: overall pressure 
change across the condensing test section as a function of flow rate; integrated 
or overall two-phase friction pressure drop multiplier as a function of satura- 
tion temperature; and local two-phase friction pressure drop multiplier as a 
function of saturation temperature. The pressure change data in each case are 
based on measurements of the two-phase potassium temperatures. The saturated 
temperature-pressure properties given in Reference 69 were used to relate the 
measured two-phase temperatures to the corresponding saturation pressures. 


Pressure Change Versus Flow Rate 

Assuming that at the condenser inlet the potassium is dry saturated vapor 
(X = 1.0), the pressure change between the inlet and a position at distance JL 
downstream from the inlet where the local vapor quality is X is given by 


P/7 - P- = 


*1 


P f g o 


where 




l 

A 

P 


r (1 ' x) 



1-x 

K 


+ 


X 


( 8 ) 


(9) 


For total condensation from X = 1 at the test section inlet to X = 0 at 
the outlet (£ = h Q ) , the corresponding pressure change, from Equation (8), is 


P - p - _T_ 
0 1 P f g o 


P . fL 

— - 1 2 0 

P 2D 1-0 


( 10 ) 
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was 


\ 


The integrated two-phase friction multiplier in Equation ( 10) , 
calculated for this portion of the data treatment using the Local, two- 
phase friction multiplier based o«n the* homogeneous ‘flow model. (K- 1) > 
which is (Reference 0) . 



1 + X 

<W - p] 


1 + X 

«W - ‘Jj 

The integration of Equation (11) to obtain 

*i-o = b /*/ :j 

x(£> J 

d t 


(ID 


( 12 ) 


is done assuming a linear variation of quality with length from the condenser 
inlet, which can be expressed simply as 



c 


The assumption of a linear variation in vapor quality is an approximation to 
the actual distribution. The comparison of calculated quality distributions 
given in Appendix E indicates that the error in local quality due to this 
approximation is a maximum of about 10% to 15% at the middle of the test 
section and about 5% at the measuring stations, compared to qualities calculated 
by energy balance. 

Overall condensing pressure change results obtained from Test Set No. 5 
(5/8-inch ID tube with instrumented 1/4-inch 0D tubular insert) are given in 
Figures 28 and 29 for vapor inlet temperatures of 1200°F and 1300°F, respectively. 
Overall pressure changes from condenser inlet to outlet calculated using 
Equation (10) are also shown in the Figures for comparison with the data. The 
analytical and experimental results are in good agreement for potassium flow 
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rates up to about 70 lbs/hr, which corresponds to vapor inlet velocities 
of 1000 ft/sec at 1200° F and 600 ft/sec at 1300°F for Test Set No. 5. The 
deviation between the experimental and calculated pressure changes becomes 
larger as the potassium flow rate increases. 


At the highest flow rates, for both 1200°F and 1300°F inlet vapor 
temperatures, the data show an overall pressure drop from condenser inlet to 
outlet, rather than the pressure increase calculated using Equation (10). 

The pressure drop across the condenser increases with increasing flow rate. 

At the lower flow rates, below about 70 Ib/hr for 1200°F (Figure 28) and 
below 100 Ib/hr for 1300°F (Figure 29), the experimental data show an overall 
pressure rise across the condenser, which is in reasonably good agreement with 
the trend calculated using Equation (10). 

Overall Two-Phase Friction Pressure Drop Multipliers 

Equation (10) can be rearranged into Equation (14), which gives the 
overall or integrated two-phase friction pressure drop multiplier 0 as a 
function of the measured overall pressure change from condenser inlet to 
outlet (Pq - Pj) and the other flow variables. 


0 


- — 
1-0 - fL 



P 

v 


1 - 


p o- p i 
o 

G /P g 

f *o J 


(14) 


Experimental overall two-phase friction pressure drop multipliers 
calculated using Equation (14) from the data of Test Sets No. 4 and 5 are 
presented in Tables 8 and 9, respectively, and are plotted in Figure 30. 
Two-phase friction pressure drop multipliers calculated using Equation (14) 
from the data of Test Set No. 1 are presented in Table 10 and are plotted 
in Figure 31. The potassium properties for these calculations are those 
given in Reference 69 evaluated at the arithmetic average temperature 
between the test section inlet and outlet. The following smooth-tube 
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friction factor equation of Blasius was used for calculation of the 
single-phase liquid friction factor (Reference 61). 

f = 0. 316/ (N (15) 

Ke 

As suggested in Reference 10, the maximum helical velocity and 
maximum helical path length were used in the calculations for the helical 
insert geometry of Test Set No. 4. 

Also shown on Figures 30 and 31 are values of the integrated two-phase 
friction multiplier calculated using the homogeneous flow model of Reference 6, 
given by Equations (11-13) and values obtained from the Martinelli model 
(Reference 31) as modified for potassium in Reference 6. Referring first to 
Figure 30, it is apparent that reasonably good agreement exists between 
calculated values using the Martinelli model and experimental values of the 
overall two-phase friction multiplier for both the helical insert data of 
Test Set No. 4 and the tubular insert data of Test Set No. 5. The experimental 
two-phase multipliers in general are higher than those calculated using the 
homogeneous model (K = 1) . 

The experimental two-phase friction multipliers obtained from the Test 
Set No. 1 data (5/8-inch ID tube without insert) are in good agreement with 
values calculated from both the Martinelli model and the homogeneous model, 
as shown in Figure 31. Nearly all of these data are in between the values 
calculated using the two analytical models. 

Local Two-Phase Friction Pressure Drop Multipliers 

Local friction pressure drop multipliers were evaluated from measured 
potassium axial temperature distributions of Test Sets No. 4 and 5. The local 
pressure gradients were determined from the measured potassium temperature 
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gradients using the potassium saturation pressure-temperature properties 
given in Reference 69. This was done for both the top and the bottom 
measuring stations in the test section. The local pressure gradients were 
determined using in each case the potassium temperatures measured with insert 
thermocouples at three positions bracketing the measuring station, and a 
least-squares procedure was used to calculate a best fit gradient through 
the three temperatures. 

The local friction pressure gradient (dP/d/) fr was estimated for each 
data point by subtracting a calculated local pressure gradient due to 
momentum change (dP/d/ ) M0M from the measured pressure gradient (dP/d^), as 
expressed by Equation (16). 

r dP dP dP, 

= dZ - ( dZ } ( 16 ) 

fr " *' MOM 

The local momentum pressure gradient can be expressed as 
f dP G 2 d l 

^’mom * 31 ( r> <17 > 


where — is given by Equation (9). 

P 

Homogeneous flow (K = 1) was assumed for the calculations, for which 
Equation (17), combined with Equation (9), reduces to 


(") 

V 


MOM 


2 

— (— 1 x dX 

p v ' pf *1 


(18) 


The local vapor quality gradient dX/djL for use in Equation (18) was 
calculated from an energy balance using the flow rate and the local heat 
flux at each measuring station. 
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The local two-phase friction pressure gradients (dP/dj£) fr , calculated 
from Equation (16), were then used in the following Equation (19) to estimate 
the corresponding local two-phase friction pressure drop multipliers 0^£ : 



2 So D P f 


fG 


tJP 

‘dZ* 


fr 


(19) 


The single-phase liquid friction factor f was calculated from Equation (15). 
For the helical insert data of Test Set No. 4, the maximum helical velocity 
was used, as suggested in Reference 10. 


Values of the local two-phase friction pressure drop multipliers 0/ 
calculated from the data of Test Sets No. 4 and 5 using Equation (19) are 
presented in Tables 8 and 9 and Figures 32 and 33 respectively. Values of 0£ 
calculated using the homogeneous flow model (K = 1) , given by Equation (11), and 
the Mart inelli flow model, modified for potassium (Reference 6), are shown on 
the Figures for comparison with the data. 


As can be seen from Figures 32 and 33 the experimental values of <t>£ 
scatter widely and agreement with calculated values using either of the two 
flow models is poor. In general, the experimental values of 0£ for the 
bottom measuring station where the vapor quality was low are higher than the 
values calculated from the two flow models, but at the top measuring station 
where the vapor quality was higher agreement with the calculated values is 
better. The vapor qualities indicated in Figures 32 and 33, X = 0. 2 at the 
bottom station and X = 0.8 at the top station, are within about 10% of the 
qualities calculated by energy balance (Appendix E) , which is sufficiently 
accurate for purposes of comparing multipliers calculated theoretically with 
those obtained from the data. 


The large scatter in the data shown in Figures 32 and 33 indicate that the 
errors in determining the local multipliers are too large for the data to be 
useful other than to show order-of-magnitude trends. 
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Overall Pressure Change (P n -P T ), psi 



Figure 28. Overall Condensing Pressure Change Versus Flow 

Rate For Test Set No. 5 At Potassium Vapor Inlet 
Temperature of 1200 F ( 5 / 8 -inch ID Tube With 
l/4-inch OD Tubular Insert) 
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Overall Pressure Change (P n -P T ), psia 



Figure 29 . Overall Condensing Pressure Change Versus Flow 

Rate for Test Set No. 5 sit Potassium Vapor Inlet 
Temperature of 1300°F ( 5 / 8 -inch ID Tube With 
l/4-inch OD Tubular Insert) 
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Figure 30. Overall Two-Phase Friction Pressure Drop Multipliers For 
Condensing From Approximately 100$ to 0$ Quality As A 
Function of Potassium Temperature Calculated From Data 
Of Test Sets No. 4 and 5 Using Equation (l4) 
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Figure J 2 . Local Two-Phase Friction Pressure Drop Multipliers As A 
Function of Potassium Temperature, Calculated From Test 
Set No. 4 Data (5/8-inch ID Tube with Helical Insert, 





VII STABILITY OF LIQUID LEVEL IN CONDENSER TUBE 


During the condensing heat transfer tests the potassium liquid level 
was always maintained in the head-tank located about one-foot downstream of 
the condenser. An alternate mode of operation is to maintain the liquid 
level in the active heat transfer zone of the condenser. An investigation of 
the fluid temperature behavior and steadiness of operation when the liquid level 
was moved from the head-tank into the condenser was done at the end of Test Sets 
No. 4 and No. 5. Results from Test Set No. 5 (5/8-inch ID tube with instrumented 
1/4-in OD tubular insert) are presented following. These results are essentially 
the same as those obtained in Test Set No. 4 (5/8-inch ID tube with instrumented 
helical insert , p/D ± = 6). 

Prior to raising the liquid level from the head-tank into the condenser 
at the end of Test Set No. 5, the following steady-state conditions were 
established: 


Boiler Power 

Average Condensing Heat Flux 
Test Section Power 
Potassium Flow Rate 
Vapor Inlet Temperature 


25 KW 

120,000 Btu/hr-ft 2 
20 KW 
65 lb/hr 
1340°F 


The test was started from these steady-state conditions by gradually reducing 
the liquid potassium pumping rate from the head-tank, while the boiler electrical 
power and the sodium flow rate were both held constant. Since the boiler power 
was held constant, the potassium vapor generation rate in the boiler and the 
potassium vapor flow rate into the condenser were nearly constant. Thus, with 
the liquid potassium flow rate having been reduced from the original steady-state 
value, the potassium vapor flow rate into the condenser became larger than the 
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liquid removal rate from the head-tank. The result of this unbalance in 
flow was to cause the liquid level to rise until it entered the condenser. 

The results of this test are presented in Figure 34, which is composed 
of two segments of a continuous oscillograph recording made during the test. 

About one-hour had elapsed from the start of the test until the beginning 
of Segment-1 in Figure 34, by which time the liquid potassium pump flow rate 
had been reduced to 56 lb/hr. Entry of the liquid level into. the condenser 
section is indicated in Segment-1 by the decrease in the potassium outlet 
temperature T from the initial 1339°F down to 1230°F over a period of several 

. K o 

seconds . 

The liquid level continued to move upward in the test section and about 
30 seconds later it reached the insert thermocouple located 5-inches upstream 
of the condenser exit. The liquid level reaching this position is indicated 
by a reduction in the potassium temperature at that point, as shown in Segment-1. 

A subsequent slight reduction in potassium liquid flow rate to 53 lb/hr 
at constant boiler power resulted in the liquid level moving gradually higher 
in the test section until it had passed the insert thermocouple located 10-inches 
upstream of the condenser exit, as shown in Segment-2. No further changes in 
the potassium liquid flow rate were made for one-hour and during this period the 
liquid level remained in the 5-inch long region of the test section between the 
insert thermocouple located 10-inches upstream of the condenser exit and the next 
higher insert thermocouple located 15-inches upstream of the exit. 

The active condensing length during the last phase of the operation, 

measured from the potassium vapor inlet end of the test section, was between 

21-inches and 26-inches long. The corresponding average condensing heat flux 

2 

and potassium vapor flow rate at steady-state were about 150,000 Btu/hr-ft and 

o 

53 lb/hr, compared to 120,000 Btu/hr-ft and 65 lb/hr, respectively, at the 
start of the test. As shown in Segment-2 of Figure 34; the insert thermocouples 
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on each side of the liquid level (10-inches and 15-inches upstream of the 
condenser exit, respectively) indicate that the local potassium temperature 
was reduced about 140°F in the 5-inch long region in which the liquid level 
was located. 

After the liquid level had been maintained in the upper position for one- 
hour, the potassium pump flow rate was increased to approximately the original 
value of 65 Ibs/hr to lower the liquid level back into the head-tank. Throughout 
the entire test only the potassium pump flow rate was varied and all other 
system control variables, such as the boiler electrical power, the sodium flow 
rate and air flow rate to the sodium cooler, were left constant. 

The following observations can be made from the test results shown in 
Figure 34: 

(1) The liquid level was maintained in a steady position in the veritical 
condenser with only small fluctuations of fluid temperature and flow rate, as 
shown in Figure 34. 

(2) A temperature differential in the potassium approximately equal to 

the difference between the potassium saturation temperature and the sodium inlet 
temperature was established over a 5-inch length of the condenser across the 
liquid level. 

(3) The liquid level could be gradually moved from one position to another 

in the condenser by adjusting the* potassium pump flow rate, without any significant 
fluctuations in temperature or flow rate other than the change in potassium axial 
temperature distribution associated with the change in potassium condensing and 
sub-cooled lengths. 
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VIII APPLICATION TO CONDENSER DESIGN 


The forced convection potassium condensing heat transfer coefficient 
data presented in Section V are applicable to design of condensers 
m which the potassium is condensed inside tubes. This application is 
discussed briefly below, primarily in reference to a 7-tube NaK-cooled, 
shell-and-tube potassium condenser. This condenser, which is described 
m more detail in Reference 71, has a tube size, shell-side geometry, 
power per tube and operating conditions which are typical of those anticipated 
for space power system condensers. 


The relationship between condensing heat transfer rate q , potassium 

condensing temperature T R , the NaK cooling fluid inlet and outlet temperatures 

"^NaK anci ^NaK ’ an< * the re< l u i re d condensing heat transfer area inside the 
I 0 

tubes A c is given by the following equation. 


q c = U i A C 


T - T 
NaK Q NaKj 


/In ( 


T - T 
K NaKj 

T - T 
K NaK 


( 20 ) 


The overall heat transfer coefficient U., based on the tube inside diameter, 
is calculated using 


1 




D. 

l 


D 

o 


h NaK 


( 21 ) 


The wall heat transfer coefficient h w in Equation (21) is given by 


1^ 

h 



w 




1 


( 22 ) 
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and the condensing heat transfer area A^ in Equation (20) for uniform 
diameter tubes is given by 


A C = ^ D i V 


(23) 


Consider the c^ign of a i i t ? W.tC v. >1-1 ooti tenser for the following 
assumed design conditions (Reference » 1; : 


Condensing heat transfer rate 
Pot assium saturate on t emperat v. 
N aK inlet t erope rat v w *e 
N aK out let t emp e t j re 
Tube out side d i ante 1 e r 
Tube wall thickness 
Tube material 

Tube spacing- 1 o--di ame* er ratio 


75 KW 
13CO°F 
1150° F 
1250° F 
0. 69-inch 
0. 035-inch 
SS 316 
1.3 


For these conditions the NaK heat t ransfer coeifi iient on the shell-side 

h , the tube wall heat transfer coefficient h and the potassium condensing 
NaK * 

heat transfer coefficient h aye estimated t o t : 

c 

h = 6, COO B t a/h :r- f t * - r F 

NaK 0 

h = 5,300 Btu/hr- ft ~ - n F 

w o 

h = 15, 000 Btu/hr- ft ' °F 

c 

The NaK heat transfer coefficient was ■;*!■ -vi ated from the equation of Dwyer 

and Tu given in Reference 40, and the wall h^at transfer coefficient was 

calculated using Equation (44). The pc tar sir,!? condensing heat transfer 

coefficient h = 15,000 Btu/hr-ft*- rt F is *:• d value based on the 

c 

data from Test Set No. 5 (5/8*-ineh It* * :.t '• with i nstrumented 1/4-inch OD 
tubular insert). As shown in Figure 2 1 thi ir.l ' e oi h^ gives a reasonably 
good approximation of the Test Set So. 5 oat a o.er the temperature range 
from 1 1 00° F to 1400° F (average deviation from the data less than + 30%). 


4 “ 


The corresponding overall heat transfer coefficient calculated using 

2 

Equation (21) is = 2470 Btu/hr-ft -°F. Using this value in Equation (20) 
together with the fluid temperatures listed above to obtain A , and sub- 
stituting A c into Equation (23) results in an estimated condensing length 
for the condenser of = 1.0 foot. Estimates of the relative sensitivity 
of the calculated condensing length to uncertainties in the heat transfer 
coefficients are as follows. 

2 

If h were actually 7500 Btu/hr-ft -°F, or 50% lower than assumed, the 

c 2 
overall heat transfer coefficient XJ would be 2120 Btu/hr-ft -°F, for which 

a condensing length of 1.17 feet would be calculated. A reduction of the 

average NaK temperature from 1200°F to 1185°F would cause the condenser to 

transfer the specified 75 KW in a 1-foot condensing length. If h were 

2 c 
22,500 Btu/hr-ft -°F, 50% higher than assumed originally, the overall heat 

2 

transfer coefficient would be 2620 Btu/hr-ft -°F and the corresponding 
calculated condensing length is 0. 95-foot. An increase in average NaK 
temperature from 1200°F to 1205°F would cause the condenser to transfer the 
75 KW in a 1-foot condensing length. This example illustrates the relative 
insensitivity of the condenser design to errors in the average condensing 
heat transfer coefficient. 

Based on the data from Test Set No. 5 (Figure 23), a reasonably 

conservative design choice for the condensing potassium heat transfer 

2 

coefficient is h^ = 10,000 Btu/hr-ft -°F. The calculated condensing length 

for the 7-tube condenser using this value is = 1.12-feet, which is 12% 

longer than the length calculated using the approximate mean value for the 

2 

Test Set No. 5 data of h = 15,000 Btu/hr-ft -°F. 

c 7 

The overall pressure change of the potassium across the condensing 
section (P Q - P^)^ can be calculated using Equation (24), derived in 
Section VI . 


- 75 - 



1 


(24) 


(P, 


Vc = 



fL 

c 

2D 


0 


1-0 


Use of the Martinelli flow model (Reference 31), modified for potassium 

(Reference 6) gives a reasonably good approximation of the overall two-phase 

friction pressure drop multiplier 0^ as can be seen by comparison with the 

data from Test Sets No. 4 and 5 shown in Figure 30 and the data from Test 

Set. No. 1 shown in Figure 31. The plots of versus potassium saturation 

temperature T for the modified Martinelli flow model in Figures 30 and 31 
K 

can be used directly to determine 0^ Q for design calculations. 
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IX CONCLUDING REMARKS 


The condensing potassium heat transfer and pressure change results 
presented provide a reasonable basis for thermal design of potassium 
condensers. The data cover the range of condensing temperatures from 
1100°F to 1400°F, local heat fluxes from 30,000 to 300,000 Btu/hr-ft 2 , 
inlet vapor Mach numbers from 0.1 to near 1.0, and local vapor qualities 
of about 0.2 and 0.8, respectively. 

Of the data presented, the condensing heat transfer coefficient data 
obtained with the instrumented helical and tubular inserts from Test Sets 
No. 4 and 5, respectively, are the most reliable, since for these tests 
the local potassium temperature was measured directly by thermocouples in the 
inserts (Figures 20-23). 


The condensing tests were conducted with the potassium in vertical 
downflow in order to eliminate the gravity force normal to the tube wall. The 
good agreement between the data taken with a straight tubular insert (Test Set 
No. 5) and those taken with a helical insert (Test Set No. 4), for which the 
radial acceleration on the fluid ranged up to 100 g f s, suggests that the effect 
of body forces such as gravity on the forced convection condensing heat transfer 
coefficient is small (Figure 22). 

All of the condensing heat transfer and pressure change data were obtained with 
the potassium liquid level maintained downstream of the condenser. Subsequently, 
the liquid level was brought up into the active condensing section to investigate 
the stability of this alternate mode of operation. 

The range of condensing temperatures tested, 1100°F to 1400°F, is 
anticipated to be the range of principal interest for design of space power 
systems using potassium as the working fluid. As discussed in Section V, there 
is some evidence that at temperatures lower than this range the condensing 
potassium heat transfer coefficients might be less than the values obtained from 
these tests (Figure 24). 
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Major conclusions indicated by the test results are as follows: 

(1) The condensing potassium heat transfer coefficients are relatively 
high, typically in excess of 10,000 Btu/hr-ft 2 -°F. The measured coefficients 
are nearly constant over the ranges of vapor saturation temperatures, mass 
velocities and heat fluxes tested (Figure 23). 

(2) The measured condensing heat transfer coefficients decrease slightly 
with increased liquid film Reynolds number (Figures 20 and 21) and with 
decreased vapor quality (Figure 23) . 

(3) Although the measured local condensing heat transfer coefficients 
are high, they are lower than coefficients calculated by considering only the 
thermal resistance due to heat conduction through the liquid film. Treatment 
of the data using an analysis based on the concept of an additional thermal 
resistance at the vapor-liquid interface (Appendix G) results in an empirical 
value of the condensation coefficient of O c = 0.2 giving a reasonable correlation 

of the data (Figure 24) . 

(4) Overall two-phase friction pressure drop multipliers calculated using 
a homogeneous flow model (equal liquid and vapor velocities) and using the 
Martinelli model, modified for potassium, are in reasonably good agreement 
with measurements using the 5/8-inch ID test section with an instrumented 
helical insert (p^ = 6) and with an instrumented J-inch OD tubular insert. 

The modified Martinelli flow model gives the best agreement with the 
experimental data for the overall two-phase friction pressure drop multiplier 
0^ (Figures 30 and 31). 

(6) Test operation included conditions for which there was a net rise 
in the potassium pressure across the condenser and conditions for which there 
was a net pressure drop (Figures 28 and 29). There was no difference in 
stability of operation between these two kinds of conditions. 
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(6) There was no observed difference in stability of operation between 
having the potassium liquid level located in the active heat transfer section 

of the test condenser and having it located downstream of the condenser in the 
head-tank (Figure 34). 
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APPENDIX A 


LIQUID HEAT TRANSFER RESULTS 


At the beginning of the experimental program, liquid heat transfer 
tests were made to gain familiarity with the facility and the test section 
instrumentation. Local liquid potassium heat transfer coefficient data 
obtained from these tests are given in Table 11 and are plotted in Figure 35. 
A detailed listing of the experimental data from which these results were 
obtained is given in Reference 4. The data cover the following range of 
vari ables: 


Reynolds Number, N 

Re 

Pec let Number, N, 

9 Pe 

Prandtl Number, N^ 

9 Pr 

Potassium Temperature, T 

K 


19,000 to 84,000 
80 to 351 
0.004 

700° F to 810°F 


The test section geometry used for these liquid coefficient tests was the 
same as was subsequently used for Test Set No.l-B (5/8-inch ID tube without 
insert), shown in Figure 7. The potassium was in vertical upflow for the liquid 
coefficient tests and was cooled by sodium in cocurrent flow on the shell-side. 


The local heat fluxes and bulk fluid temperatures at the measuring stations 
were calculated from energy balances using the measured flow rates and fluid 
inlet and outlet temperatures. The overall heat transfer coefficient for the 
test section was assumed to be constant for the energy balance analysis. The 
error introduced by this assumption was small since more than half of the thermal 
resistance between the two fluids was due to the thick-walled nickel tube. The 
inner and outer tube wall temperatures were determined by a least-squares 
procedure to fit the calculated heat fluxes and the measured temperatures inside 
the tube wall at the measuring stations to the integrated Fourier heat conduction 
equation, as described in Appendix E. 
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The local liquid potassium heat transfer coefficient h R at each 
measuring station was then calculated from the inner wall heat flux q/A^ 
the inner wall surface temperature T and the local potassium bulk fluid 
temperature T using Equation (A-l). 

A 



Potassium fluid property values given in Reference 66 were used for these 
calculations . 

Equations (A-2) , (A-3) and (A-4) are plotted in Figure 35 for comparison 
with the data. 

N = 5.0 + 0.025 N 0,8 (A-2) 

Nu Pe 

N = 0.625 N °' 4 (A-3) 

Nu Pe 

N = 3.66 + 0.0055 N„ (A-4) 

Nu Pe 

The experimental results are all below values calculated using the theoretical 
Equation (A-2) of Seban and Shimazaki (Reference 50). The experimental values 
are generally less than those calculated using the empirical Equation (A-3) 
of Lubarsky and Kaufman (Reference 47) , but the general trend with Peclet number 
appears to be about the same as that given by Equation (A-3). The author’s 
empirical Equation (A-4), which is similar in form to Equation (A-2), provides 
a reasonably good fit to the data for Peclet numbers above about 150. At 
lower Peclet numbers the data fall below values calculated from Equation (A-4) . 

The data obtained at the bottom measuring station (>L/ D i = 19) are 
believed to be more reliable than those obtained at the top measuring station 
(JL/ D a = 38), for the following reasons: 
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(1) The difference between the potassium bulk temperature and the 
wall temperature was greater at the bottom measuring station, 
thereby making the coefficient determinations at the bottom 
station. less sensitive to thermocouple error. 

(2) Four tube-wall thermocouples were functioning at the bottom 
measuring station during the liquid coefficient tests, whereas 
only three thermocouples were functioning at the top measuring 
station, which allowed a more reliable determination of the inner 
wall temperature at the bottom measuring station to be made. 

Thus, due to these uncertainties the small dependence on /// D. indicated 
by the data in Figure 35 is not necessarily true. 

For runs 1-28 (Table 11 ) a maximum error in the experimental Nusselt 

numbers of + 15% is estimated. A slightly larger error is possible for 

data runs 29-31, due to the fact that at the low Peclet numbers for these 

runs (Np^ < 120) the difference between the measured sodium and potassium 

temperatures was small, less than 2°F at the test section outlet. For the 

rest of the data (N_ -^120) the smallest difference between the measured 

Pe 

fluid temperatures at the test section outlet was 4.5°F. 

The sodium Reynolds number was maintained in the range from 21,000 to 
29,000. No trend in the potassium Nusselt number as a function of sodium 
flow rate is discernible in the data, which indicates that the dependence of 
the sodium heat transfer coefficient on Reynolds number was small over this 
range. An average sodium Nusselt number of 5.3 was obtained for the sodium 
Peclet number range of 113 to 156. The measured data varied from this average 
value with a standard deviation of 31%. Agreement between the average experimental 
value of the sodium Nusselt number of 5.3 and the corresponding value of 5.6 
calculated using the equation of Dwyer and Tu for an annulus (Reference 63) is 
good. 
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Nusselt Number 



Curve-a: Equation A-2, Seban and Shimazaki (Reference 50) 

Cu'rve-b: Equation A-J, Lubarsky and Kaufman (Reference 47) 

Curve-c: Equation A -4, author's correlation 


Figure 35* Liquid Potassium Heat Transfer Data Obtained 
in 5/8-inch ID Tube at 700°F to 8lO°F. 
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APPENDIX B 


ELECTROMAGNETIC FLOWMETER CALIBRATION 


Electromagnetic flowmeters were used in both the sodium loop and the 
potassium loop. Both flowmeters were used for the liquid heat transfer 
coefficient tests described in Appendix A. For the condensing tests, 
however, the potassium loop flowmeter was not used, due to the small potassium 
flow rates in condensing, and the potassium flow rate was determined instead 
by energy balance across the test section using the measured sodium flow rate. 


The equation for calculation of the flow rate of an electrical conducting 
fluid from the emf generated by an electromagnetic flowmeter, from Reference 67, 
is: 


!° 4 D . ^ 

W/p f E ' = 3.18 FK X K 2 K 3 


(B-l) 


where 


2 

r ^ 2i 

D. 

1 + (_i) + R’ 

!- (^) 

o 

o 


The constant K 2 in Equation (B-l) is a correction for end-effects and is 
equal to 0.989 (Reference 67). The factor Kg is a function of the magnet 
temperature and dimensions and for these experiments was between 0.99 and 1.0 
The pertinent dimensions and magnetic flow ratings of the two electromagnetic 
flowmeters at room temperature are as follows: 


Flowmeter 

D. , inch 

D , inch 

—o^—— 

F, gauss 

Potassium 

0.618 

0.841 

1989 

Sodium 

0.610 

0.843 

2019 
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At the beginning of the liquid heat transfer coefficient tests described 
in Appendix A a discrepancy in heat balance across the test section between 
the potassium and sodium sides was found. The discrepancy was that the 
indicated heat transfer rate from the potassium was about 15% higher than the 
indicated heat transfer rate to the sodium. It was therefore decided to 
calibrate the potassium electromagnetic flowmeter in-place in the loop and to 
then calibrate the sodium flowmeter relative to the potassium flowmeter. 

A calorimetric method was chosen for the calibration. The thermal 
calorimeter used is illustrated in Figure 36. An immersion heater was used 
for heat input to the calorimeter to assure that all of the heat input to the 
calorimeter went into potassium sensible heat. To determine the potassium 
flow rate, three parameters have to be known: 


(1) Heat input to the calorimeter, q, Btu/sec 

(2) Specific heat of liquid potassium, c , Btu/lb-°F 

P K 

(3) The potassium temperature increase across the calorimeter 

due to the heat input, (T - T ), °F 

K 0 K I 

The heat input to the calorimeter was measured by a standard single-phase 
wattmeter having full scale ranges of 500, 1000, and 2000 watts and an accuracy 
of + 1/4% of full scale. The calorimeter heat loss and the systematic tem- 
perature measurement error between the inlet and outlet potassium thermocouples 
were determined by a technique similar to that used in the sodium thermocouple 
calibrations, described in Appendix C. 


With zero heat input to the calorimeter, the potassium loses heat to the 
surroundings as a function of the temperatures of the calorimeter chamber and 
the surroundings. For the condition of zero heat input to the calorimeter 
(case-1) the heat loss from the calorimeter is given by: 


q 


LI 


(Wc ) 

P K1 




(B-3) 


- 86 - 


where A is the systematic temperature measurement error between the inlet and 
outlet potassium thermocouples. With heat input q^ to the calorimeter (case-2) 
an energy balance gives: 


ca 


“ <1 


L2 


= (Wc ) 
P 


(T„ 


- T„ 


“ A). 


K2 


(B-4) 


Measurements with and without heat input were made at the same average 

potassium temperature, for which the calorimeter chamber temperatures and 

corresponding heat losses q L were thus also the same. With this condition 

imposed q ca can be calculated from the following Equation (B-5), which is 

obtained by combining Equations (B-3) and (B-4) with q = a 

LI m L2 * 


ca 


(T„ 


" T., 


< Wc d >K1 


(V %>K2 


(B-5) 


The potassium electromagnetic flowmeter was used to adjust the flow during 
the calibrations so that (Wc p ) R1 = (WO^. With this condition, together 
with the assumption that the systematic temperature measurement error between 
the inlet and outlet potassium thermocouples A was independent of the calorimeter 
heat input (A ± = A g ) , Equation (B-5) reduces to 


w _ ca 

K2 c ~ TIT - T ) +7t — T ) 1 (B-6) 

P K2 [ K q \ 2 ' K1 WlJ 

Equation (B-6) was employed to determine the potassium flow rate using 

measured temperatures and heat input to the thermal calorimeter. For assumed 

errors of + 0.5% in the calorimeter heat input q^, + 2.2% in the potassium 

specific heat c pR and + 1.6% in the measured potassium temperature differences 

<* K “ T k ^ 2 and (T k " T k ^ i ’ the estimated standard error in the calculation 

0 1 10 
of flow rate using Equation (B-6) is + 3.1%. 
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Flow rates calculated from Equation (B-6) were used as the basis for 
calibration of the potassium electromagnetic flowmeter over the following 
ranges of calorimeter conditions: 


Potassium flow rate 
Potassium Temperature Difference 
Average Potassium Temperature 
Calorimeter heat input 


0.14 to 0.31 lbs/sec 
11 . 5°F to 43 . 4°F 
600°F to 700° F 
0.76 to 1.6 Btu/sec 


Four sets of calibrations resulted in an average ratio of the calorimeter 
to electromagnetic flowmeter flow rates of 1.145 with a standard error of 
+ 1.3%. Multiplication of the flow rate indicated by the electromagnetic 
flowmeter by the factor 1.145 brought the test section liquid- liquid heat 
balances between the potassium and sodium sides into agreement within a standard 
error of 4 - 6.7%. 


The test section heat balance error was calculated from Equation (B 7). 


Error = 1 


q L + 


Na 


(B-7) 


The test section heat loss q was determined during the calibration of the 

L 

sodium thermocouples, as discussed in Appendix C. The heat transfer rate from 
the potassium and the net heat transfer rate to the sodium were calculated 
from fluid temperature and flow rate measurements made during liquid operation 
using Equations (B-8) and (B-9) , respectively. 


q 


K 


W c 
K pKM 




(B-8) 


q 


Na 


W c ( T 

Na pNaM Na r 



(B-9) 
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The fluid temperatures used in Equations (B-8) and (B-9) were corrected values 

based on the thermocouple calibrations discussed in Appendix C. The potassium 

flow rates W R were corrected values based on the thermal calorimeter calibration 

results discussed above. The sodium flow rates used were uncorrected values 

obtained directly from the electromagnetic flowmeter output together with 

Equation (B-l). Reference 68 was used for the mean values c „ of the sodium 

pM 

and potassium specific heats at test section temperatures. 


The standard error in the test section heat balance calculated using 
Equation (B-7) was + 6.7%. The heat balance error includes the error in sodium 
flow rate determination and errors from other sources. The heat balance error 
being small indicates that the sodium flow rate error was also small. On this 
basis, the sodium flow rates obtained from the electromagnetic flowmeter output 
together with Equation (B-l) were used directly, without further correction, 
for the condensing data reductions. The potassium flow rates for the condensing 
tests were determined by energy balance across the test section using the 
measured sodium flow rates (Appendix E) . 


- 89 - 






APPENDIX C 


THERMOCOUPLE CALIBRATION 

Stainless steel sheathed, capped, 0.040-inch diameter, chrome 1-alumel 
thermocouples were used, due to their stability and relatively high emf output 
in the 1100°F to 1400°F temperature range. The high output of the chromel- 
alumel thermocouples, approximately 23 microvolts /°F, minimized temperature 
errors due to instrumentation noise, which was generally less than + 3 micro- 
volts. 

Since the condensing heat transfer coefficients for potassium are high, 
m the order of 10,000 Btu/hr-ft 2 -°F, accurate calibrations of the test section 
thermocouples were necessary. For calculation of the condensing heat transfer 
coefficient the difference between the vapor and wall temperatures is needed. 
Therefore, in-place calibrations of all test section thermocouples relative to 
each other were done. 

Initial calibrations indicated inhomogeneities in the chromel-alumel 
sheathed thermocouples, which might have been caused by bending the thermo- 
couples during installation. These inhomogeneities were detected by moving a 
soldering iron at about 700°F along the exposed thermocouple length after 
installation and observing the change in thermocouple output. All thermocouples 
which indicated greater than 50 microvolts (2°F) output change during the 
homogeneity check were replaced. In addition an individual ice junction was 
provided for each thermocouple to reduce errors caused by junction emf 1 s 

For intercalibration of the sodium inlet and outlet well thermocouples, 
sodium fluid temperature measurements were obtained at one temperature and 
two different sodium flow rates, with the potassium loop evacuated and with 
no heating of the test section (Reference 4). From these data the test 
section heat loss and the individual thermocouple corrections were determined 
independently as functions of average sodium temperature, as shown below. For 
two runs under the above conditions at two different flow rates, Equations C-l 
and 02 are applicable: 
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(C-l) 


’Ll * <* V K.1 lT »a : - T »a 0 * A) ! 


’l 2 * <” V »-2 "»*, - T Na 0 * " 2 


where A is the required correction between 
thermocouples. 


(02) 

the sodium inlet and outlet 


At the same average sodium temperature the average shell temperatures and, 
consequently, the test section heat losses will be the same for both runs, i.e., 
q = q . Also, if each thermocouple has a systematic correction which is 

“j^l ' 

independent of flow rate, the systematic thermocouple correction between the sodium 
inlet and outlet thermocouples can be obtained from Equation A-3. 


A = 


% - T ». 0 >1 - " N ,, ~ T K. 0 >2 


l - 


<■ ' V Hal 

<" V K.2 


+ (T 


Na„ 


T ) 

Na^l 


(C-3) 


The assumption of thermocouple corrections independent of flow rate is 
reasonable for the test section geometry used. 


This method of calibration yields only relative calibrations between the 
inlet and outlet sodium well thermocouples. The correction, A, was applied to 
thermocouple readings used in the calculation of the amount of heat transferred 
to the sodium from the potassium during condensing tests. The test section 
heat loss was simultaneously determined during these thermocouple calibration 
runs. 


In a similar manner , an attempt was made to determine relative calibrations 
between the potassium inlet and outlet well thermocouples. However, these 
calibrations were unsuccessful, probably due to the lack of fluid mixing at the 
test section inlet and outlet on the potassium side. Therefore, to obtain 
relative calibrations between the potassium fluid thermocouples and the nickel 
tube wall thermocouples, a different procedure was used. 
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For calibration of the potassium thermocouples, the test section annulus 

was evacuated and potassium vapor was passed through the test section while 

the test section shell was maintained at the potassium vapor temperature by use 
of electrical heating. The small amount of vapor passing through the test 

section from the boiler was condensed by heat losses in the head-tank below the 

condenser. Test section thermocouple readings were obtained under these 
conditions for each test geometry at 1200°F and 1300»F, both before and after 
the heat transfer tests. From these data a systematic correction for each test 
section thermocouple relative to a selected standard thermocouple was obtained 
as a linear function of temperature. For the last two test sets (Test Sets No. 4 
and 5), the maximum thermocouple correction relative to the standard thermocouple 
was less than 3°F. Agreement between the thermocouple corrections determined 
before and after condensing operation was generally better than 1°F. 
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APPENDIX D 


TEMPERATURE FIELD DISTORTION DUE TO THE THERMOCOUPLE 
HOLES IN THE THICK-WALL NICKEL CONDENSER TUBE 


An analogue using teledeltos paper was used to determine the magnitude 
of the temperature field distortion caused by the axial thermocouple holes in 
the thick walled nickel tube of the test section. A flux plot was made on a 
six-times enlarged cross-section of the test section geometry. After the 
teledeltos paper was cut to size, it was checked for uniformity of resistance 
in the vicinity of the thermocouple holes before the holes were cut. The 
maximum deviation of the measured potential from the calculated values was 
+ 0.3%, with an applied potential difference of 10 volts between the inner and 
outer radii. 

The flux plot determined from the analogue is shown in Figure 37. Maximum, 
minimum and average potentials measured around the circumference of each hole 
are listed in Table 12 together with the values calculated for the undistorted 
geometry at the radial position of each hole center. The calculated and 
measured values agree within a maximum error of 3.5%. On the basis of this 
small error, the temperatures indicated by the wall thermocouples were assumed 
for the data reductions to be equal to the temperature at the radial position 
of the respective thermocouple hole centers. 

The noncircular holes shown in Figure 37 simulate those used in the first 
test section built, which was subsequently discarded. All the other test 
sections had drilled circular thermocouple holes. This slight difference in 
shape of the thermocouple holes between what is shown in Figure 37 and those 
used for the condensing tests would not have enough effect on the temperature 
field to alter the principal conclusion of this study, that the effect of the 
thermocouple holes on the measured wall temperature distributions is negligible. 
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Figure 37 • 


Flux Plot Showing, Analogue of Temperature Field In 
Thick-Walled Nickel Condenser Tube. 


APPENDIX E 


CONDENSING DATA REDUCTION PROCEDURES 


Fluid Properties 

The potassium thermodynamic properties used for reduction of the 
condensing data were obtained from Reference 69. Specific heats for 
liquid sodium and liquid potassium were taken from Reference 68. Other 
sodium and potassium liquid transport properties used for the data 
reductions were obtained from Reference 66. 


Potassium Flow Rate 

In general, the potassium flow rates used for the condensing tests 
were too small to measure accurately with the electromagnetic flowmeter. 
For this reason the potassium flow rates for the condensing tests were 
determined indirectly by calculated energy balance across the test section 
using the measured sodium flow rate and measured sodium inlet and outlet 
temperatures, as follows. 


Assuming that the potassium entered the test section as saturated vapor 
and that the effect of subcooling of the liquid film at the test section 
outlet was negligible, the potassium flow rate W was calculated from 

K 

VJ _ _ 

K ~ X (E-l) 

By an energy balance across the test section the heat transfer rate from the 
potassium q R was calculated from 

q K = q Na + q L (E-2) 

Combining Equations (E-l) and (E-2) gives for the potassium flow rate 

W K = X (q Na + q L ) (E-3) 
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I 


The test section heat loss rate q was determined by test, as described 

L 

in Appendix C. The net heat transfer rate to the sodium q was calculated 
from 


Na 


= W XT C PxT 
Na Na 


(T 


Na„ 


- T. 


Na, 


(E-4) 


The sodium temperature difference in Equation (E-4) was determined using 
measured bulk fluid inlet and outlet temperatures, corrected for relative 
thermocouple errors as described in Appendix C. The sodium flow rate was 
calculated from the sodium loop electromagnetic flowmeter output using 
Equation (B-l), as discussed in Appendix B. 


Local Wall Temperature and Heat Flux 

The local temperatures at the inner surface of the condenser tube wall T 
and the local heat fluxes q" were determined from temperature measurements using 
thermocouples located inside the wall at the top and bottom measuring stations, 
respectively. The arrangement of the wall thermocouples is described in 
Section III. 


The Fourier heat conduction equation (Reference 12), Equation (E-5), 


R i R 

T = T . - qV — In (|-) 
w wi H i - R. 

k 1 

w 


(E-5) 


was fitted to the measured internal wall temperatures T by a least-squares 

w 

fit procedure for each data run. The local temperature at the inner surface 

of the wall was then calculated from the fitted Equation (E-5) , as illustrated 

for two typical cases in Figure 16. The calculations were done using the 

following matrix equation, Equation (E-6) , derived from Equation (E-5) to 

express the inner wall surface temperature T in terms of the least-squares 

Wi 

fit to the measured internal wall temperature T^ at the respective radial 
positions R . 
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(E-7) 


Inside wall temperatures T wj were normally measured at five different 


locations 


inside the wall at both measuring stations, for which the summation limit in 


Eauations (E-6) and (E-7) is N =5. 

o 


The local heat flux at the inner wall surface was calculated for each case 

using Equation (E-8) , derived from Equation (E-5) to express the local inner 

wall surface heat flux ql* in terms of the least-squares fit to the measured 

internal wall temperatures T at R , 

wj j’ 


k 

M w 

q i = r?9 
1 


N 


N 


r 


In (v±) 


N 

r 


i 

N 


£ 


wj 


R i 

T . In (^1) 
wj V R. 

l 


(E-8) 
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in which 0 is given by Equation (E-7) . 


The mean value of the tube wall thermal conductivity k^ used for these 
calculations was evaluated for each case at the mean wall temperature T , 
which was calculated using 



The thermal conductivities of the INCO Nickel-270 condenser tube used 
were values measured by Battelle Memorial Institute, Columbus, Ohio, given 
in Figure 38. 

Local Potassium Saturation Temperature 

The local potassium temperature axial distributions were measured 
directly using the insert thermocouples for Test Sets No. 4 and 5. The local 
potassium temperature at the measuring stations for Test Sets No. 1, 2 and 
3, for which no insert thermocouples were used, was determined by linear 
interpolation between bulk potassium temperature measurements at the test 
section inlet and outlet. An examination of the dependence of the experimental 
condensing heat transfer coefficients on the potassium temperature distribution 
is given in Appendix H. 

For Test Set No. 4 an instrumented helical insert (p/D^ = 6) was used. 

For the data from Test Set No. 4 the insert thermocouple temperatures were 
corrected by analysis to account for the radial pressure rise and corresponding 
saturation temperature rise in swirl flow to obtain the estimated saturation 
temperature at the tube wall. The analysis used is as follows. 


In cylindrical co-ordinates the equation of motion is 


V 

dp g. 

dR g R 
o 


(E-10) 
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I 


With the assumption that the axial component of the velocity , V is 
independent of the radius, Equation (E-ll) is obtained. 


dP 

dR 



2 n V 



2 


R 


(E-ll) 


If the liquid fraction in the two-phase mixture is neglected, and the 

density is assumed to be the vapor density, Equation (E-ll) can be integrated 

to obtain Equation (E-12), which gives the difference between the pressure 

at the condensing surface, P , and the pressure at the insert centerbody P . 

D DCB 


P - P _ ( 2L, 2 (JL) 2 

D DCB V 2g ' P V P/D' 
o v 



(E-12) 


The temperature measured by the helical insert thermocouples was assumed 

to be the saturation temperature corresponding to the static pressure at the 

insert centerbody, P . The radial pressure difference due to the swirling 

DCB 

flow, as calculated from Equation (E-12) was then added to P to give P , 

DCB D 

the pressure existing at the liquid-vapor interface. A corrected value of the 

potassium saturation temperature T corresponding to P was then determined and 

K D 

used to calculate the condensing heat transfer coefficient. 


The effect of this correction was to decrease the experimental condensing 
heat transfer coefficients for the helical insert data of Test Set No. 4, as 
a result of increasing the vapor saturation temperature from the measured 
centerbody saturation temperature to that corresponding to the pressure at 
the liquid-vapor interface. With this procedure, potassium saturation 
temperature corrections up to 20°F were obtained. As indicated in Figure 22, 
reasonably good agreement between the data for the 5/8-inch ID tube with helical 
insert (Test Set No. 4) and the data for the tubular insert (Test Set No. 5) 
was obtained after application of the correction. 
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Condensing Heat Transfer Coeffi ci ent 


Experimental values of the local condensing heat transfer coefficient 
h^ were calculated for each measuring station using Equation (E-13). 



T ,) 

wi 


(E-13) 


and the values of T ^ and q" obtained from Equations (E-6) and (E-8) . The 
corresponding Nusselt condensing ratio ^ UQ and the liquid film Reynolds 
number N Ref used in the correlations in Section V were calculated from 
Equations (E-14) and (E-15), respectively. 
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(E-14) 


4W ( 1-X) 

N = - — 

Ref tt D i U 


(E-15) 


Fluid properties used in Equations (E-14) and (E-15) were those of the saturated 
liquid (Reference 66). 


Local Vapor Quality 

For Test Set No. 5 (5/8-inch ID tube with instrumented 1/4-inch OD tubular 
insert) the local vapor quality at the measuring stations, listed in Table 3 and 
used in Equation (E-15) , were calculated by energy balance using the measured 
potassium temperature distribution, inlet and exit sodium temperature, sodium 
flow rate, and the experimental condensing heat transfer coefficients. The 
procedure used for this calculation was to first estimate a value for the overall 
heat transfer coefficient using calculated values of the shell-side sodium heat 
transfer coefficient and tube wall heat transfer coefficient together with the 
measured condensing heat transfer coefficient. Then, the heat flux and sodium 
temperature distributions were calculated using the estimated overall heat 
transfer coefficient, the measured potassium temperature distribution and the 
sodium exit temperature. The calculation proceeds step-wise from sodium exit to 
inlet. The difference between the calculated and measured value of the sodium 
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inlet temperature is used as the closure criterion. If agreement within about 
^ F is not obtained on the first iteration the value of the overall heat transfer 
coefficient is adjusted and the procedure repeated. This process is continued 
until the measured and calculated sodium inlet temperature agree within about J°F. 
Finally, the vapor quality distribution was calculated by energy balance on the 
fluid using step-wise integration of the heat flux proceeding from the potassium 
vapor inlet. 

The results of the calculation for one of the data runs is shown in 
Figure 39. As can be seen from Figure 39 the difference between the vapor 
quality calculated by energy balance and the quality estimated by linear 
interpolation is about 5% at both measuring stations. The measured potassium 
temperature distribution for the data shown in Figure 39 (Run No. 37), is one 
of the most non-linear obtained from Test Set No. 5 (Figure 17). Thus, the 
differences between the calculated quality and quality estimated by linear 
interpolation shown in Figure 39 should be among the largest for Test Set No. 5. 

A similar procedure for calculating the local quality was attempted for 
the Test Set No. 4 data (5/8-inch ID tube with instrumented helical insert, 
p/D i = 6), using an additional step in the analysis to correct for the effects 
of the swirl-flow on the saturation temperature at the wall. The results of 
the calculation for one of the data runs is shown in Figure 40. For this run 
the difference between the vapor quality calculated by energy balance and the 
quality estimated by linear interpolation is about 5% at the bottom measuring 
station and about 8% at the top measuring station. Since the required correction 
for swirl flow effects adds additional uncertainty to the local vapor quality 
calculations and since the differences between calculated qualities and qualities 
estimated by linear interpolation (Figure 40) are small, linear interpolation was 
used for estimating the local quality at the measuring stations for use in 
Equation (E-15) for the Test Set No. 4 data. 

For the data of Test Sets No. 1, 2 and 3, the local quality at the 
measuring stations was estimated for use in Equation (E-15) by linear interpo- 
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lation. Based on the calculated quality comparisons shown in Figures 39 
and 40 the estimated error in local vapor quality at the measuring stations 
is about 5%. The error in estimating the local vapor quality causes a small 
error in the calculation of the local liquid film Reynolds number. The error 
in the vapor quality estimates is not relevant in determination of the 
experimental condensing heat transfer coefficients. 

It was assumed for the estimates of local vapor quality that the potassium 
enters the condenser as saturated vapor and exits as saturated liquid. The 
errors introduced by these approximations are assumed to be negligible. The 
piping from the boiler to the condenser was well-insulated to prevent any 
significant heat losses from the vapor. The potassium exit piping was 
similarly well-insulated. Estimates of the heat losses between the condenser 
exit and the liquid level in the head-tank downstream of the condenser (Figure 1) 
indicate that the vapor quality of the potassium at the condenser exit was less 
than 1%. The vapor inlet piping from the boiler to the condenser was well 
insulated to maintain the inlet vapor quality at about 100%. 

Pressure Change 

The pressure change data presented in Section VI were obtained from 
measurements of the two-phase potassium temperatures. The potassium pressures 
were assumed to be the saturation pressures corresponding to the measured local 
fluid temperatures. The fluid property data in Reference 69 were used for the 
pressure change data reductions. 

Liquid Film Coefficient For Calculation Of Vapor Phase Resistance 

The liquid-film heat transfer coefficient h f , defined by Equation (5), was 
calculated using Dukler* s model (Reference 28) for use in Equation (6) to obtain 
the vapor-phase resistance 1/h^ from the Test Set No. 5 data. The Nusselt 
condensing ratio, (h^/k^) (V^/g^) , calculated using Efcikler’ s model is given 

in Figure 41 as a function of liquid-film Reynolds Number and T*, the 
dimensionless shear stress at the vapor-liquid interface. 
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The dimensionless shear stress T* is calculated using 


T* = 
V 


p ( <vX > 1/3 


(E-16) 


for which the vapor shear stress is given by 


dP. 


T v = » (D i - V <dZ> v 


(E-17) 


The local pressure gradient of the vapor due to friction was estimated using 


, _ f p V 

,dP. _ v K v v 

V v ' 2« 0 <■>, - V 


(E-18) 


where 


V 

v 


4 "k x 


n p. 




(E-19) 


The vapor-phase friction factor was calculated from the Blasius smooth-tube 
equation (Reference 61), 


f 

v 


0.316 


(n y 

Rev 


(E-20) 


Although Dukler 1 s thermal analysis of the condensing problem is in error 
for low Prandtl number fluids, as has been shown by Lee (Reference 27), the 
Dukler model accounts for shear at the vapor-liquid interface and appears to 
correlate liquid-film thickness data reasonably well (References 33, 34 and 35). 
For this reason Dukler f s film thickness analysis was used for the calculations. 


-105- 



o 



- 106 - 


Figure 38 . Thermal Conductivity of INCO Nickel 270 Used for 
Condenser Tube as Measured by BMI. 




Local Quality, X, dimensionless 



Figure 39. Calculated Potassium Quality Distribution Compared 
to Linear Interpolation for Run No. 37 of Test Set 
No. 5 (5/8-inch ID Tube with Instrumented l/4-inch 
OD Tubular Insert) 



Local Quality, X, dimensionless 


\ /Station 


■ ✓Linear 

Interpolation 


Calculated by 
Energy Balance^ 


Run No. 21 
Test Set No. 4 


Bottom 

Station 

V 

V \ 

\\ 

\\ 

\\ 


10 15 20 25 30 
Distance from Test Section Inlet, inches 


35 36 


Figure 40. Calculated Potassium Quality Distribution Compared 
to Linear Interpolation for Run No. 21 of Test Set 
No. 4 (5/8-inch ID Tube with Instrumented Helical 
Insert, p/D, = 6) 
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APPENDIX F 


HEAT TRANSFER COEFFICIENT ERROR ANALYSIS 

The error in the experimental local condensing heat transfer coefficients 
is a function of the errors in measurement of the following local quantities: 

(1) Potassium saturation temperature 

(2) Heat flux 

(3) Tube inner wall temperature 

An analysis of these measurement errors was done for the data of Test Set 
No. 5 (5/8-inch ID tube with instrumented 1/4-inch OD tubular insert). 

T emperature Error - Two sources of temperature measurement error were 

considered, random measurement errors and errors due to thermocouple drift 
alone. 


To determine the probable random error in the measurement of the potassium 
saturation temperature, readings from three of the insert thermocouples were 
considered. Three repetitive readings from each of the three thermocouples for 
three different groups of runs were used for the statistical analysis. The 
thermocouple identification numbers (T/C No.), the distance of the thermocouple 
from the condenser inlet (/, inches) and the probable errors calculated from 
the thermocouple readings are listed below. 


Insert 
T/C No . 

23 

28 

29 




inches 


1 

26 

31 


Probable Error 
+ 0. 1°F 
+ 0. 4°F 
+ 0. 4°F 
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The identification numbers for the three repetitive runs in each data 
group, as listed in Reference 9, and the dates on which each data group 
used for these calculations were taken are: Runs 19, 20, 21, taken 2/27/65; 

Runs 25, 26, 27, taken 2/28/65; and Runs 28, 29, 30, taken 3/2/65. The 
average of the probable errors for these three thermocouples is + 0.3°F, 
which was assumed to be the probable error in the potassium thermocouple 
temperatures E,^. 

The probable random errors for two of the nickel tube wall thermocouples, 
one at the top measuring station and one at the bottom measuring station, were 
also estimated from readings taken during these same runs. The thermocouple 
identification numbers (T/C No. ), measuring station at which the thermocouples 
were located, the radial position of the thermocouples in the tube wall 
(R, inches) and the probable errors calculated from the thermocouple readings 
are listed below. 


Wall 
T/C No. 

Station 

R, inches 

Probable Error 

11 

Top 

0.409 

+ 0. 4°F 

16 

Bot t om 

0.408 

+ 0. 3°F 


The average of the probable errors for these two thermocouples is + 0.3°F, 

which was assumed to be the probable error in the wall thermocouple temperatures 

E m • 

Tw 


The thermocouple temperatures used were corrected values, based on the 
calibrations discussed in Appendix C. To make allowance for thermocouple 
drift with time, the thermocouple corrections used were the averages between 
the calibrations performed at the beginning of test operation and those 
performed at the end of test operation. The differences in indicated temperatures 
from the thermocouples between the before-test calibrations and the end-of-test 
calibrations were assumed to be drift. The absolute values of the drifts 
listed below were measured at 1300°F. The errors due to thermocouple drift were 
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assumed to be one-half the total drift between the before-test calibrations 
and the end-of-test calibrations. 

Thermocouple Number 



11 

16 

23 

28 

29 

Drift, °F 

+ 1.6 

+2.4 

+ 1.0 

+0.6 

+0.4 

Error, °F 

+0.8 

+ 1.2 

+0.5 

+0.3 

+0.2 


The probable error due to drift for the five thermocouples is assumed to be 
the average of the five individual errors. On this basis the probable error 
due to thermocouple drift E^ was calculated to be 0.6°F. 


Heat Flux Error - The error in the heat flux q" is due to uncertainties in the 

tube wall thermal conductivity and uncertainty in the measured temperature 

gradient in the condenser tube wall (Appendix E) . For this analysis these 

two error sources were assumed to have probable errors of 2% each, which results 

2 

in an estimated probable error for the heat flux of E „ = 0.028 q" , Btu/hr-ft . 

q i 


Wall Temperature Error - The probable error in the temperature at the tube 
wall inner surface T was assumed for this analysis to be 1% of the total 
temperature drop across the nickel tube wall. Based on an estimated average 

2 

heat transfer coefficient for conduction through the tube wall of 1300 Btu/hr-ft -°F, 
relative to the tube inner diameter (5/8-inch ID), the estimated error in the 
tube wall inner surface temperature is E Twi = 0.01 q"/1300, °F. 


Condensing Coefficient Error - The probable error in the condensing heat transfer 
coefficient was estimated from these component errors using Equation (F-l), 

2 / / \ 2 2 i 


<-*-> + 
- «■; 




T, 2 2 

E + E 
TK d 

T - T . 

K wi 


h £ 


2 E 2 
TW + d 


T - T 
K wi 


E, 


Twi 


T - T . 
K wi 


(F-l) 
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for which, in summary, the component probable errors used were: 


(1) 

V 

thermocouple drift error = 0.6°F 

(2) 

e tk> 

potassium thermocouple error = 0.3°F 

(3) 

E Tw’ 

wall thermocouple error = 0.3°F 

(4) 

E Twi * 

tube surface temperature error = 0.01 q M /1300°F 

2 ^ 

(5) 

E", 

n 

heat flux error = 0.028 q" , Btu/hr-ft 


q 


At typical test conditions of q" = 150,000 Btu/hr-ft 2 and (T - T ) = 15°F 

1 K wi 2 

which corresponds to a condensing heat transfer coefficient h^ = 15,000 Btu/hr-ft -°F, 

a probable error in h c of + 15% was calculated using Equation (F-l) . Assuming 

a random distribution of errors, this value corresponds to a standard error in 

h of + 22% (Reference 70). 
c - 

Figure 42 shows values of the probable error in the condensing heat transfer 
coefficient as a function of heat flux and pot assium-t o-wall temperature difference 
estimated using Equation (F-l). As shown in the Figure, the probable error 
decreases with increased temperature difference at constant heat flux, and 
increases with increased heat flux at constant temperature difference. 

A similar error analysis was not performed for the data taken in the other 

Test Sets. It is believed that the data from Test Sets No. 1, 2 and 3, for which 

direct measurements of the local potassium temperature at the measuring stations 

were not made, have a larger error than the Test Set No. 5 data, due to additional 

errors involved in estimating the local potassium saturation temperature T , as 

K 

discussed in Appendix H. The data from Test Set No. 4 (5/8-inch ID tube with 
instrumented helical insert, p/D^ = 6), for which the local potassium temperature 
was measured directly, should have approximately the same errors as estimated 
in this Appendix for Test Set No. 5 plus an additional uncertainty due to the 
correction for effects of swirl flow on the saturation temperature at the tube 
wall, as discussed in Appendix E. 
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Probable Error In Condensing Heat Transfer Coefficient 


W lx: 



Figure 42. Estimated Probable Error In Condensing Heat 
Transfer Coefficient As Function Of Fluid- 
To-Wall Temperature Difference And Heat Flux 
For Test Set No. 5 Data ( 5 / 8 -inch ID Tube With 
Instrumented l/4-inch OD Tubular Insert) 



APPENDIX G 


ANALYSIS OF POTASSIUM VAPOR-PHASE THERMAL RESISTANCE 


In order to compare with theory the values of the vapor-phase coefficient 
°k^ a: *- ne d from the Test No. 5 data using Equation (6), an analysis was 
conducted based on a procedure suggested by the work of Schrage (Reference 32) 
and used by Rohsenow and Sukhatme (Reference 16). The analysis is as follows. 


From kinetic gas theory (Reference 36) , the mass flux of molecules 
impinging on a surface is given by 

Mg 

G m = P (— — ) (G-l) 

2TTRT 

When equilibrium exists between a saturated vapor and its liquid, the 
mass flux of molecules leaving the liquid surface equals the mass flux of vapor 
molecules condensing on the liquid surface. That is, the net rate of heat 
transfer is zero. This equality is expressed by 


P 

v 



2ttR T 


v 



i 

) 


2ttR T 

S 


(G-2) 


in which: 


? v = pressure of the saturated vapor 
T^ = temperature of the saturated vapor 

T g = temperature of the liquid at the liquid-vapor interface 

P = saturation pressure corresponding to T 
s s 

When a net rate of mass transfer occurs, such as in condensing, a condition of 
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I 

! 

nonequi librium exists at the vapor-liquid interface. For condensing, the 

net mass flux of molecules G leaving the vapor and condensing into the liquid 

c 

can be expressed as 


G 

c 


a P 

C V 



h 

) 


2 TT R T 

v 


a p 

e s 



i 

) 


2 n R T 

s 


(G-3) 


where 

a = condensation coefficient, the mass flux of molecules from the 
c 

vapor which actually condense divided by the mass flux predicted 
by Equation (G-l) for P v and 


a 

e 


evaporation coefficient, the mass flux of molecules which actually 

leave the liquid surface divided by the mass flux predicted by 

Equation (G-l) for P and T . 

s s 


It was assumed for the treatment of the Test Set No. 5 data that 
With this assumption Equation (G-3) reduces to Equation (G-4) . 



Substituting from 


P 

v 




2 n R T 


v 

Equation (G-l) 


P 

s 




2TTR T 

s 


into Equation (G-4) gives 


G 

c 



G 

ms 


(G-4) 


(G-5) 


Assuming that all of the heat transferred from the vapor to the liquid 
film is by condensation, an energy balance on a unit area of the heat transfer 
surface, using Equation (G-5) , gives 
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(G-6) 


q 


ft 


= a 


c 




X g - X 

v mv s 


G 

ms 




= a A (^ G ) 
c m 


where the difference A is between the two sets of pressure and temperature 

conditions P , T and P . T . 

v v s s 

A vapor phase heat transfer coefficient h can be defined as 

v 

h v = q " /AT v (G-7) 

where q" is the condensing heat flux and AT y is the corresponding difference 

between the bulk vapor temperature T and the liquid surface temperature T . 

v s 

Substituting into Equation (G-7) from Equation (G-6) for q" results in 


h = a 
v c 


A (^ G ) 
m 

AT 

v 


(G-8) 


For small values of the temperature difference AT , the latent heat of 
vaporization can be assumed to be constant and the righthand side of Equation (G-8) 
can be expressed as a derivative to give Equation (G-9) 


dG 

h = a X — E 
v c dT 

v 


(G-9) 


Substitution from Equation (G-l) into Equation (G-9) results in Equation (G-10) 
for the vapor-phase heat transfer coefficient. 


h = o X 

V c 



2 n R T 


(^ - ^) 
v dT 2T 


(G-10) 
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Equation (G-10) was used to calculate the vapor-phase heat transfer coefficients 
for comparison with the values obtained from the Test Set No. 5 data (Section V, 
Equation 7). Saturation properties for potassium from Reference 69 were used for 
these calculations. 

Values of the vapor-phase coefficient for potassium as a function of saturation 

temperature with the condensation coefficient & as a parameter, calculated using 

Equation (G-10), are shown in Figure 43. As can be seen from the Figure, the 

calculated vapor phase heat transfer coefficient increases with increased 

saturation temperature T and is proportional to the condensation coefficient G . 

K 

It is interesting to compare the relative importance of the vapor-phase 

heat transfer coefficient for potassium with that of a higher Prandtl number fluid 

such as water. Figure 44 shows the vapor-phase heat transfer coefficients 

calculated for water and potassium, respectively, using Equation (G-10) with 

a = 1.0. The corresponding liquid-film heat transfer coefficients calculated 
c 

using Nusselts model (Reference 12) for an assumed liquid-film Reynolds number 
3 

N =10 are shown in the Figure for comparison. 

Ref 

Figure 44 shows that over the pressure range from 0.1 to 100 psia, the vapor 

2 

phase heat transfer coefficient for water at ® = 1.0 is in the order of 10 to 

4 ^ 

more than 10 times as large as the liquid film heat transfer coefficient. A 

value of a =0.3 has recently been reported for water (Reference 37). Thus, 
c 

the vapor phase thermal resistance for water is very small relative to the thermal 

resistance due to conduction through the liquid film, even for condensation 

coefficients as small as a =0.1. In contrast, for potassium, even with 

c 

condensation coefficients as large as ^ = 1 . 0 , the vapor phase coefficient is 

about equal to the liquid film coefficient at 1 psia (T gat = 990°F) and is less 
than ten times the liquid film coefficient at 25 psia (T ga ^ = 1500°F) . 
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Calculated Vapor Phase Heat Transfer Coefficient, h , Btu/hr-ft 
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Figure 


Potassium Saturation Temperature, T K , °F 

Vapor Phase Condensing Heat Transfer Coefficient Calculated 
For Potassium From Kinetic Theory of Gases Using Equation (7) 
For Various Values of 6 . 



Heat Transfer Coefficients for the Vapor Phase, h , and the Liquid Film 
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Comparison of Liquid Film Coefficients Calculated by Nusselt's 
Model (Reference 12) with Vapor Phase Coefficients Calculated 
by Equation (7) for Condensing Water and Condensing Potassium. 
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APPENDIX H 


DEPENDENCE OF EXPERIMENTAL CONDENSING COEFFICIENTS 
ON POTASSIUM TEMPERATURE DISTRIBUTION 


As indicated in Appendix E, linear interpolation between the bulk 
potassium temperatures measured at the condenser inlet and outlet, respectively, 
was used to estimate the local potassium temperature at the measuring stations 
for the data of Test Sets No. 1, 2 and 3, The effect on the heat transfer 
coefficient accuracy of the linear potassium temperature distribution approxima- 
tion can be studied by comparing coefficients obtained from Test Set No. 5 
(5/8-inch ID tube with instrumented 1/4-in. OD tubular insert), which are based 
on direct measurement of the local potassium temperature, with values calculated 
from the same data using linear interpolation to estimate the local potassium 
temperatures . 

Figure 45-a shows the condensing heat transfer coefficient data from Test 
Set No. 5 obtained using the measured local potassium temperatures. these are 
the same data as presented in Figure 21 except that the three low-temperature 
data points (discussed in Section V) have been deleted for this comparison. 

The same data, but with the coefficients obtained using local potassium temperatures 
estimated by linear interpolation from measured condenser inlet and outlet 
bulk potassium temperatures, are shown for comparison in Figure 45-b* 

Comparison of Figure 45-b with Figure 45-a indicates that using linear 
interpolation to estimate the local potassium temperature causes increased scatter 
of the heat transfer coefficients compared to the coefficients obtained using the 
measured local potassium temperatures. The increased scatter is particularly 
evident for the data taken at the top measuring station ( ^/D. = 10). 
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TABLE ^ 

CONDENSING HEAT TRANSFER DATA FROM TEST SET NO. 5 
(5/8-inch ID Tube with Instrumented 1/4 -inch OD Tubular Insert 
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TABLE 4 


VAPOR PHASE COEFFICIENTS CALCULATED FROM TEST SET NO. S DATA 
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CONDENSING HEAT TRANSFER DATA FROM TEST SET NO 
(5/8-inoh ID Tube wi'.'nout Insert) 


U 


t'-rnvo ov- 4 - 4- mmo nvo 


me— c— b-ovovovovovoo 

CMCMCMCMCMCMCMCMCunm 


00000000000 
C-VO 1 — l t-COvp mcO O ovp 
mCJvOO ONO OVVO CM O rA 


vo tpi^r in^r cr\ 4 -vo 00 e— 0 
voChCMO OVD 4 e-t—OCO 

in n cr\ n in^t mvo vo cm n 
mvo vo O-t-e- t— e-r— COCO 




ao 

cr 3 ] 

n 

■p 


ml 

o 

a 


X 

*H | 

a> 

ao 

\ 

n 

h 



O til 

sz ^ 


00000 
ino HHin 
in CM OVCMOO 


COCO CM 04 " 


4- 4 - O CJV o 


nnnn m 

00000 


O Q OQ O 
vo in cm 04-. 
mm4 o t- 




00000 000 

c— 4 -vo 4 - in rHco o 
mm rH e-co rH ovoo 

4 -nD-mcM e-004- 

e-t— mitco coco o\ 


mov 00004 - m^r m 

cm vo vo o o m«H ov 

in^t CM CM H O 0 4 
HHHHCM CM CM CM 


co h e-mm 4 o o\ 
o ovco m cm mmm 
m 4 m n m m 4 m 
00000 000 


00000 000 

00 C-if WHO b-4 

rH m 00 m4 00 cm in 

Oo' t KN CM* rH rHinCM^ 


vo 04 h kv rncvi in 

, it 4 VO lAt- 4 CM H 
5*4 H 1 — I rH CM CM CMCMrn 

Eh O I HrlrlrlH HrlH 


(X, 

O 

g 


a 

X 


rH m t—4 cm in mao o m4- 

OHHHOOCMHCM 4 CVI 

mnmmmmmnmcMCM 

00 000000000 


31 


00000000000 

t-b-CMOt-e-OC-CM 4 CO 

t— rH n cm t-aovo m4 e-ov 


t-nao movvo rnrnao cm m 
X pt, me- f-ao e- O' — t' — — c — 00 co 

Eh O I HHrlrlHHrlrlHH rH 
rlHrlrlHHHHHHH 


rno tAHio inoo4-cocM 
m lAt-Koo c-e-t— c- e- e-co 

E *4 PH rH rH.rH HHHHHHHH 
£hO rlHHiHHrlrlrtHHH 


OOOOO 
CM O mco CM 
t-CM mm 4 

00 c^c^oooo 


lAH rH CM rH 

VO 4 e-VO O 
CM CM CM 4 LH 
CM CM CM CM CM 


CTVVO CTv OV 4 
44 H ova 0 
CM CM in CM CM 
OOOOO 


ovovo mao 
e— 4 e— 4 0 




a) 

af 5 


OOOOO 000 
e— ovcm rncM t-o in 
rH CM - 4 “ CM - 4 " 4 int- 

e^tnmoo rH <— T rn 
inm- 4 - - 4 - e— t-e— co 


4 m 4 mt\ 4 iacm 

no CM CM O CM CM vo 
ON CO P— CX\ C 7 \ t— 4 GO 
CM CM CM CM if if it in 


mavmmmt— cm4 


00000 000 


OOOOO 000 
vo <-h r - cm m onn 
avao cm cm cm ovme- 

O VO CM nit rH 00 o 
CMrHrHrHrH H H (\j 


OiTiNlAOiAinOOHH 
rH if O n n- 4 " H 4 H H H 
00 00 CTv On O O h rH CM nit 
OOOOHHHHHHH 


nnnnnnnnnnn 

VOVOVOVOVOVOVOVOVOVOVO 


mmmmmmmmmmm 


X Pel 

£h o I 


rH rH OVVO 00 
WHO H H 




X Pel 
Eh o I 


00 rAO OVrH OVO CM 

n mvo 4^4 hh 

H rH rH. CM CM CM CM m 




tH O 1 


CM rH 4 OVOV CO rH o 


nCMVO CM OV OVOvO 

n n cm cm n nit it 


mocono avit rn 
nnmit t- nn rn 

H H rH CM CM CM CM n 
H H H H H rHrIrt 


c • 

3 0 O rH 

K s rn cm nit mvo e-co a\ >-» <-h 


Q> OOOOO 

e o cm mcM n 

tH nO rH cm CM 

Eh CM O O O O 


mo -4- mm o mo 
it o nit it o n n 
ovocMoOav rH cm vo 

rH CM O O O rH rH rH 


nnnnn 
vovo vo vo vo 


■4* it it it 4“ it 4 4" 1 

vo vo vo VO VO vo vo^ 


e-e-e-e-e- e-ne— 


c • 

3 O 

oc 55 rn cm nit m 


c • 

3 o 

pc zl 


rn cm nit m vo t— ao 


- 129 - 



vo| 

a 

9 

EH 


Q 

W 

CQ 


EH 

< 

W 

ffl 

o 

s 

H 

CO 

S 

B 

S 

O 

o 


o o o 

VO mo 
a aco 


o o o o 
o o o o 

CO AGO A- 

00 A ACO 
OH^-KV 

I- 1 n I r— I r— I 


o o 
o o 
cvi m 


o o o o o 
o o o o o 

a a Aptoo 
vo tnKvt^rA 

t— A- H O A 
r — l I — l OJ CVI C\J 


oooooooo 

OOOOOHQOiA 
A-PT CM vo ApJ- C-CM 

vo a- aco cm acapt 
Aih in cm ocovovo 

CM CM H H H 


O CA A- 

co a-pt 
O KNO 

CM CM CM 


aco a rH i 

LAO A- A 
O CA ACC i 

A CM A A A A pj--=j- avovo VO A A A CM rH rH r 


O A AVO ov 
f- AGO O A- 


pj-vo CAO OVO A CD 
CM «H ApT AA-pTVQ 
vo VO OG I— i st Ovvo PT 




Q 


CM O VO 

CM VO A 00 
A 00 A- CM 

A- rH 

caca ah o 





iH 

o A 

ST pf A 00 o 



N 


AO O 

CM O CArH 

o ca 

A AO 00 O 



z 

rH O O 

H CM H rH 

rH O 

rH O <H O rH 



C 

z 1 

1 1 






o 

p 






-p 

-p 

o, 






u 

a} 

-P 

CM 1 
-P 






0 

CO 

CO 

Or 

o o o 

O O o o 

o O 

O O O o o 



1 

COCO rH 

A- CAvO A- 

Pf A 

rH vo A- 00 CA 


6 

u 

VO Pf VO 

AHH^ 

pf o 

pf A 00 CM A 

CVI 

£ 

o 

O JC 





M 

4-> 

sz \ 

VO CM rH 

Apf O CA 


a- crvvo opt 
pf rH A A A 


-P 

3 

Pf 1 i 

4 AAA 



o 

-P 





o 

G 

CQ 

CQ 






•rH 








Ph 



CM VO VO 

rH A- CO A- 

rH A 

A CAVO ArH 

EH 



W Ph| 

O ACM 

O O O CA 

CM ON 

[—CACAO O 



£h O 1 

1 CM CM A 

A CM CM CM 

pf a 

CM rH CM Pf Pf 

W 

TJ 



i — 1 i — 1 rH 

rH rH rH H 

rH pH 

H i — 1 rl rl H 

CO 

0 








G 


CM 

-P 





Eh 

0 






C0 

Ph 


<H 

OO o 

o o o o 

o o 

o o o o o 



1 

aovooo 

CM rH o O 

AO 

o o o o o 

H 


u 

rH C — rH 

VO 00 rH rH 

vo pj- 

a- o t — pt vo 

Eh 

Eh 


:_l £ 









rH f — t~ — 

A A-C0 A- 

o CM 

VO o O Apf 
Apt a-vooo 




3 

VO VO oo 

cooooo 

O A 

s 

T5 


-P 

PQ 


*H rH 

rH rH 

i — 1 rH rH ( — 1 rH 

o 

0 






pc: 

P 

c8 







G 



pf pf o 

pf pf A O 

00 CM 

ovo 

0 

o 

a> 

CM rH O 

A- CM VO 00 

00* OJ 

00 o AAA 

< 

B 

y 

A- AH 

vo o cm av 

rH pf 

^vo AAA 
COCO O rH CM 
i — 1 i — 1 rH 

§ 

3 


z K 

! ApT A- 

l 

A AVO VO 

a-ca 


4 AO A-A-VO O CM 
pf rHVO-=rO lAAA 
CAO AVO CA A O CM 
O 1 I rH rH CM CM (A lA 


o o 
oo 

CO rH 

vo 00 

00 A- 
CM CM 


VO A 
pj- <j\ 
CM CA 


p 1 -vo 

O A 
00 rH 
O rH 


OOOOOOOO 
[ — rH CM A-VO CAO VO 
A-VO CA rH 00 Pt A- rH 

CM A A - 00 CM O 00 A- 
A Apt AO CAO CO 


A CT\A- Apt H rH VO 
st vo APf A A A rH 
ACM CM CM CM CM CM CM 
H i — I i — IrHi — I < — I i — I H 


O O 
CM 00 
AVO 


A- A 
rH A 

■=r -=r 


oooooooo 

OOOVOAAOVA 
ACACA A-OO A CO VO 


CM CM CAVO VO H 

VO H H C\J O 04 CM 
CM CM 00 VO Pf VO O A- 
CM O vo ApT A ACM 


A- (A 
CM 00 
apt 


P 

o*' 








w 

rH 


A ACM 

o A-avoo 

rH OV 

sT CD sT CM Ov 

rH 00 A CO St AAA 

GO CM 



00 A A- 

vooo rH A- 

OVO 

CAVO A- CM H 

VO 00 CA CA CA A-VO CM 

rH O 

G 

ii 


CTnCM rH 

OCO A A 

ACM 

A 00 VO CM A 

VO H CM PT VO Pf t 

O CA 

H 



o o o 

sr O CM A 

CM ^ 

rH O rH AO 

OrHiHCM-HCMPTCM 

VO rH 

t 

Q 

S5 

1 1 1 

1 l l l 


1 1 rH 

A 1 1 1 t 1 1 1 

CM 

C 







1 

i 

o 









525 

C 

Ph 








o 

o 







s: 

-p 

CM 

oo o 

o o o o 

O O 

O O O O O 

OOOOOOOO 

o o 

p 

flj 


oo cm 

O A- ACM 

ACM 

00 rH rH ACM 

C0 CA CA 0J 00 00 A A 

ovo 


-p 

Cm 

. CTnCM o\ 

O A A A 

a a\ 

O CAvO A- A- 

VO A-VO A CAVO r i A- 

Pf vo 


co 

1 







5 


o u 

St 00 A 

CM rH CM 00 

Avo 

VO o 00 O rH 

st rH A00 CA A-OOvO 

r-IPf 


Or 

jc x: 

A 

sT A00 00 

00 St 

A A A rH vo 

vo -=r -4- 00 A00 ACA 

CAvO 


o 

\ 

1 1 1 

rH 1 1 rH 

rH 

1 1 rH A 

O 1 1 1 1 1 rH 1 

co 

0 

Eh 

3 

P> 


1 1 




1 

rQ 


m 





1 


3 

r ■ 









c3 



. 00 CM A 

CJN rH CM A 

O A 

00 CA ACM o 

CM A CVI CA rH 00 CAp(- 

A A 




ovovoo 

O', O O ON 

CM CT\ 

VU Omj\ o o 

pfVOAAACMCVJH 

rH A 

Q 


W pH 

rH CM A 

CM CM CM CM 

a 

CM rH CM -4- Pi- 

ACMCMCMCMCVlCMCM 

pr PT 

H 


Eh O ! 

I rH rH rH 

H rl H rt 

rH rH 

rn rH rH rH H 

rHi — 1 i — It — 1 rl rl H i — 1 

rH rH 

s: 









o 


u 

A-O A 

rH CO OV A 

CM CO 

CM vo CM rH 

OVO A-Pf rH A- OCO 

CM Pf 

G 



rH A A 

A- CD rH CM 

A- CO 

O CACA A^-O 

H A4 rH 00 00 CM 00 

ACA 

•rH 



A 

A A C — A- 

VOCD 

CACO O O CM 

CM O A-VO Pj- A ACM 

PT A 



jQ 




i — 1 i — 1 rH 

rH rH 

rH rH 

CO 









\ 









m 



00 CM 00 

CT\ O O A 

O A 

A 00 CM OJ O 

CM AOCO OC0C0A 

CM CM* 



H 

1 OvCT\ CT\ 

CvO O CTv 

CM CA 

VO CA CA O O 

pTVDAAACMCMH 

rH A 



Ph 

rH CM A 

CM CM CM CM 

-4" A 

OJ rH CM P)- 

ACM CM CM CM CM CM CM 

-=r ■=* 



Eh O 

1 rH rH rH 

H rH rl rl 

i — l i — l 

H rl rl rl H 

rlrlrlrlHrlrlrl 

r— 1 rH 


OO o 
oo o 

OAOV 


AO O A 
4 O AH 
A- A CO CM 
O |H rH CM 


AO 
p}- A 
CM CA 
O rH 


O O O AO 
CM O AH O 

8 A rH A- CM 
o rH rH CM 


O o o A AO O o 
CM ApT 4 H Oj CM O 
CM O AVO 00 CA O rH 
O rH rH rH rH H CM CM 


O o 

AO 

Pf rH 


A- 


■=r 

oo 


rH cm A pj- avo A- coca 


O rH CM Apf 


AVO C— 00 CA O rH. CM 
r — 1 i — I H H H 0J CM CM 


ApT 
CM CM. 


130 


TABLE 7 

CONDENSING HEAT TRANSFER DATA FROM TEST SET NO. 
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Data calculation not valid for Runs No. 2, 7 8 iii and 

due to subcooled liquid at potassium outlet ’ ' 



TABLE 8 


TWO-PHASE FRICTION PRESSURE DROP MULTIPLIERS OBTAINED FROM TEST SET NO. 4 DATA 
(5/8-inch ID Tube With Instrumented Helical Insert, p/D^ ” 6) 


Run. , . 

T 

KI 

T 

KO 

W K 

v i 

"TV 

0 (2) 
1-0 

^(3) 

^(3) 


°F 

°F 

lbs/hr. 

ft/sec. 

psi 


Top Station 

Bottom Station 

1 

1182.8 

1179.1 

23.3 

380 

. 105 

1319 

1854 

873 

2 

1247.1 

1245.1 

21.0 

234 

.078 

954 

1379 

639 

3 

1155.6 

1145.9 

40.6 

783 

.231 

1622 

2124 

1268 

4 

1221.8 

1216.1 

38.8 

501 

.208 

1156 

1625 

1026 

5 

1295.8 

1292.4 

35.8 

306 

. 168 

823 

1359 

835 

6 

1399.6 

1398.4 

29.6 

150 

.091 

510 

559 

346 

7 

1391.7 

1390.4 

46.7 

246 

.099 

466 

543 

347 

8 

1282.4 

1278.6 

53.4 

490 

.177 

811 

1025 

662 

9 

1207.7 

1198.2 

57.9 

812 

.300 

1258 

1580 

1012 

10 

1109.9 

1067.2 

54.5 

1418 

.705 

2462 

3544 

2018 

11 

1146.9 

1128.8 

50.0 

1020 

.401 

1832 

2552 

1461 

12 

1203.0 

1195.4 

47.5 

684 

.234 

1256 

1682 

1007 

13 

1211.3 

1200.9 

57.0 

783 

.334 

1259 

1654 

1011 

14 

1207.0 

1182.3 

61.8 

870 

.746 

1518 

2153 

1135 

15 

1260.0 

1235.5 

69.8 

725 

.967 

1275 

1807 

895 

16 

1295.8 

1283.3 

69.4 

593 

.618 

965 

1297 

694 

17 

1407.6 

1405.1 

61.5 

300 

.198 

485 

607 

384 

18 

1405.4 

1401. 1 

83.1 

410 

.341 

520 

625 

431 

19 

1301.0 

1269.6 

87.3 

725 

1.49 

1164 

1674 

799 

20 

1255.5 

1199.9 

87.6 

933 

2.00 

1522 

2303 

1082 

21 

1212.5 

1058.0 

89.9 

1225 

3.45 

2196 

3697 

2144 

22 

1267.0 

1140.1 

107.5 

1074 

4.07 

1778 

2889 

1535 

23 

1306.5 

1249.0 

104 

839 

2.63 

1317 

2066 

931 

24 

1405 . 2 

1380.5 

107.5 

531 

1.87 

837 

1205 

545 

25 

1401.3 

1370.0 

121 

608 

2.30 

859 

1216 

556 

26 

1310.3 

1218.4 

121 

957 

3.96 

1425 

2275 

1111 

27 

1321.0 

1166.9 

142 

1063 

6.08 

1529 

2358 

1393 

28 

1407 

1364.9 

131.5 

644 

3.10 

923 

1328 

623 

29 

1395.5 

1323.7 

144.6 

747 

4.74 

1102 

1620 

725 

30 

1142.4 

1130.9 

34.8 

731 

.253 

1807 

2570 

1443 

31 

1110.4 

1091.2 

37.7 

978 

.362 

2272 

3370 

1836 

32 

1093.2 

1081.8 

21.4 

624 

. 187 

2426 

3684 

1694 

33 

1152.5 

1147.0 

18.8 

369 

. 130 

1736 

2665 

1150 


(1) Operating conditions for each run are given in Table 2 

(2) 0 calculated from data using Equation (14) 

(3) calculated from data using Equation (19) 
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TABLE 9 


TWO-PHASE FRICTION PRESSURE DROP MULTIPLIERS OBTAINED FROM TEST SET NO. 5 DATA 
(5/8-inch ID Tube with Instrumented 1/4-inch OD Tubular Insert) 



T 

KI 

°F 

T 

KO 

°F 

lbs/hr. 

v i 

ft/sec. 

<W Vo< 2 > 

psi 

% 1 2 (3) 

Top Station 

0 / (3) 

Bottom Station 

l 

1106.2 

1109. 1 

23.0 

585 

.053 

2419 

2821 

2262 

2 

1143.1 

1144.8 

22.7 

453 

.039 

1896 

2118 

1742 

3 

1194.7 

1195.2 

20.3 

294 

.015 

1391 

1539 

1315 

4 

1257.4 

1257.4 

18.9 

190 

0 

1025 

1231 

945 

5 

1252.7 

1253.3 

' 34.5 

356 

.033 

1125 

1608 

1494 

6 

1200.8 

1203.6 

36.1 

503 

.088 

1433 

1778 

1764 

7 

1142.3 

1151.2 

38.4 

710 

.206 

1939 

2108 

2044 

8 

1112.6 

1122.0 

38.1 

928 

. 184 

2513 

2979 

2516 

9 

1112.7 

1128.5 

44.7 

1089 

.316 

2494 

3207 

2578 

10 

1149.4 

1160.9 

43.9 

841 

.280 

1871 

2228 

2024 

11 

1197.2 

1204.2 

44.2 

630 

.219 

1376 

1548 

1500 

12 

1301.3 

1303.4 

42.7 

338 

. 109 

763 

755 

795 

13 

1404.4 

1404 . 9 

41.0 

194 

.040 

492 

515 

397 

14 

1402,4 

1402.8 

45.8 

219 

.032 

544 

622 

461 

15 

1307. 1 

1309.9 

51.1 

392 

. 145 

783 

837 

820 

16 

1238.8 

1245.5 

53.0 

592 

.259 

1108 

1180 

1221 

17 

1206. 1 

1215.4 

53.9 

728 

.305 

1374 

1472 

1506 

18 

1186.8 

1205.9 

69.1 

1048 

.582 

1627 

1814 

1691 

19 

1250.0 

1260.0 

69.8 

731 

.410 

1124 

1169 

1201 

20 

1295.2 

1301.0 

69.6 

569 

.289 

899 

878 

954 

21 

1393.8 

1395.1 

65.0 

323 

.099 

587 

609 

574 

22 

1400.4 

1401.3 

82.3 

396 

.071 

656 

741 

561 

23 

1302 * 4 

1308. 1 

87.0 

685 

.295 

997 

1278 

834 

24 

1296.0 

1301.5 

86.9 

707 

.276 

1049 

1258 

902 

25 

1256.2 

1267.2 

87.3 

883 

.465 

1238 

1840 

959 

26 

1211.3 

1209,2 

86.0 

1126 

-.069 

1999 

4154 

1191 

27 

1260.4 

1274.4 

105.0 

1038 

.611 

1286 

1708 

1139 

28 

1312.0 

1311.8 

107.2 

802 

-.011 

1190 

2626 

3395 

29 

1413.2 

1413". 7 

102.7 

466 

.041 

680 

906 

646 

30 

1410.5 

1416.0 

121.3 

557 

.454 

588 

935 

362 

31 

1309,8 

1320.8 

123.1 

931 

. 598 

1036 

2233 

572 

32 

13120 

1308.0 

138.0 

1011 

-.429 

1338 

3790 

179 

33 

1296. 5 

1278,7 

131.0 

1064 

-.853 

1589 

4149 

263 

34 

1264 . 0 

1119.1 

126.6 

1226 

-.438 

2950 

5946 

2000 

35 

1198.0 

1074 .. 3 

108.3 

1535 

-2.75 

3369 

6431 

3007 

36 

1167.0 

884.4 

96.7 

1657 

-3.48 

4380 

10,747 

6703 

37 

1162.3 

1182.0 

68.8 

1214 

.527 

2024 

1081 

1778 


(1) Operating Conditions for each run are given in Tab: 

(2) calculated from data using Equation (14) 

(3) 0 calculated from data using Equation (19) 
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TABLE 10 


TWO-PHASE FRICTION PRESSURE DROP MULTIPLIERS 
OBTAINED FROM TEST SET NO. 1 DATA 
(5/8-inch ID Tube Without Insert) 



T 

T 


V T 

(P -P T ) 

0, A 

z: <» 

KI 

KO 

K 

I 

0 I 

1-0 

°F 

°F 

Ibs/hr 

ft/sec 

psi 


Group A 
1 

1153.4 

1161.0 

25.7 

410 

. 176 

1657 

2 

1170.0 

1176.8 

27.1 

395 

. 176 

1517 

3 

1175.3 

1181.7 

27.6 

385 

. 169 

1560 

4 

1180.7 

1186.5 

27.8 

380 

. 158 

1601 

5 

1175.5 

1182.6 

29.4 

410 

. 188 

1618 

6 

1173.2 

1180.5 

29.4 

415 

. 191 

1635 

7 

1169.4 

1177.3 

29.5 

430 

.202 

1627 

8 

1169.7 

1177.4 

29.5 

430 

. 199 

1646 

9 

1174.5 

1181.7 

29.0 

405 

. 189 

1572 

10 

1178.5 

1185.4 

30.3 

415 

. 187 

1690 

11 

1181.7 

1188.4 

30.6 

410 

. 183 

1690 

Group C 

1 

1135.0 

1148.5 

33.2 

600 

.290 

2306 

2 

1130.2 

1143.6 

32. 1 

595 

.284 

2330 

3 

1158. 1 

1165.8 

26.3 

410 

. 187 

1518 

4 

1248.9 

1250.6 

22.9 

205 

.068 

1052 

5 

1269.6 

1274.6 

39.7 

315 

.214 

894 

6 

1238.6 

1244.9 

40 

375 

.234 

1190 

7 

1214.3 

1223.2 

40.2 

440 

.292 

1271 

8 

1311.4 

1315.7 

46.0 

295 

.228 

734 


(1) Operating conditions for each run are given in Table 5 

(2) Q calculated from data using Equation (14) 


TABLE 11 


LIQUID POTASSIUM HEAT TRANSFER DATA* 




Bot t om 

Station. ^£/D, 

= 19 

Top 

Station. 

= 38 

Run 

Pe 

(T -T . 
K wi 

) 4 

q . . 10 Btu 

\u 


1« 

4 

q . , 10 Btu 

Nu 

No, 


ll 

hr-ft^ 


°F 

1 hr- ft 2 


1 

351 

12.7 

2.87 

5.09 

6.85 

1.77 

5. 79 

2 

340 

12.2 

2.76 

5.10 

6.49 

1.68 

5.81 

3 

315 

11.9 

2.65 

4.97 

6. 18 

1.57 

5.68 

4 

289 

12.0 

2.60 

4.86 

5.98 

1.50 

5.58 

5 

264 

10.3 

2.22 

4.79 

5.02 

1.22 

5.42 

6 

244 

10.4 

2.19 

4.69 

4.87 

1.17 

5.33 

7 

219 

10.7 

2. 19 

4. 57 

4.75 

1.10 

5. 15 

8 

188 

10,8 

2. 13 

4.39 

4.41 

0.99 

5.00 

9 

171 

11.0 

2. 11 

4.28 

4. 17 

0.92 

4.90 

10 

143 

10.8 

2.04 

4. 19 

3.51 

0.78 

4.93 

11 

162 

9.8 

1.97 

4.49 

3.64 

0.82 

4.99 

12 

176 

10.0 

1.98 

4.43 

4.00 

0.89 

4 . 93 

13 

202 

9.8 

1.96 

4.48 

4.30 

0.96 

4.99 

14 

233 

9.8 

1.97 

4.53 

4.73 

1.05 

4.98 

15 

256 

9 , 4 

1.99 

4.77 

4.57 

1. 10 

5.41 

16 

278 

9.3 

1.98 

4. 78 

3. 19 

1. 14 

5.34 

17 

308 

9.5 

2.11 

5.04 

5.01 

1.26 

5.66 

18 

331 

9.4 

2. 15 

5. 17 

5. 10 

1.32 

5.83 

19 

345 

9.4 

2.16 

5.24 

5. 16 

1.34 

5.90 

20 

341 

14.4 

3.13 

5 . 00 

8. 18 

1.97 

5.49 

21 

324 

14.4 

3.05 

4 c 84 

8. 14 

1.89 

5.30 

22 

302 

14. 6 

3.02 

4 , 72 

8.06 

1.83 

5. 17 

23 

270 

15. 1 

3.04 

4 o 59 

7.82 

1.76 

5. 11 

24 

239 

15.2 

2.86 

4.25 

7.86 

1.58 

4.51 

25 

238 

14.9 

2.95 

4.45 

7. 18 

1.61 

5.00 

26 

279 

14.4 

2.99 

4.67 

6.63 

1.74 

6. 11 

27 

305 

14.6 

3.00 

4. 65 

7.38 

1.82 

5. 60 

28 

336 

14.4 

3.19 

5. 02 

7.94 

1.97 

5.59 

29 

115 

6 . 6 

1.12 

3.82 

2.29 

0.37 

3.62 

30 

94 

6 . 1 

1.00 

3.67 

1.76 

0.26 

3.43 

31 

80 

6.4 

0.94 

3.33 

1.73 

0.21 

2.75 


-^Obtained with test section geometry of Test Set No. 1-B (Figure 7) , 5/8-inch ID 

thick-wall nickel tube, potassium in vertical up-flow cooled by sodium in cocurrent 
flow. 
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TABLE 12 


RESULTS FROM FLUX-PLOT ANALOGUE OF TEMPERATURE 
FIELD IN THICK-WALL NICKEL CONDENSER TUBE U) 


T/C 

Hole No. 

Radius, 

inch 

Calc' d. . 
Potential^' 

Potential 

Maximum 

From Flux Plot 
Minimum 

Analogue** 

Average 

R i 

0.313 

1000 

1000 

1000 

1000 

1 

0.398 

758 

875 

630 

753 

2 

0.449 

637 

748 

515 

632 

3 

0.518 

494 

600 

400 

500 

4 

0.601 

346 

435 

260 

348 

5 

0.698 

196 

275 

130 

203 


(1) Potentials are normalized to a value of 1,000 at the inner surface 
(radius R^) corresponding to 10-volts in the analogue test. 

(2) Calculation based on tube wall without thermocouple holes. 

(3) See Figure 36. 
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